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I. 


INTRODUCTION  AND  SUMMARY 


This  quarterly  report  describes  progress  made  during  the  past 
quarter,  1  July  1971  to  30  September  1971,  on  the  Extended  Evaluation  of  ALPA, 
NORSAR  and  LPE  Data  program  being  conducted  by  Texas  Instruments  Incorpor¬ 
ated  at  the  Seismic  Array  Analysis  Center.  The  program  consists  of  the  follow¬ 
ing  seven  tasks: 

•  Continued  evaluation  of  the  Alaskan  Long  Period  Array  (ALPA) 

•  Evaluation  of  the  long  period  Norwegian  Seismic  Array  (NORSAR) 

•  Evaluation  of  the  short  period  Norwegian  Seismic  Array 

•  Evaluation  of  the  stations  of  the  Long  Period  Experiment  (LPE) 

network 

•  Investigation  of  network  processing  and  analysis  techniques 

•  Adaptive  processing  studies 

•  Investigation  of  high- resolution  frequency- wave  number  spectral 
estimation  techniques 

The  software  required  to  perform  the  evaluation  was  developed 
under  a  previous  contract  (  Contract  F33657-60-C-1063).  During  the  past  quart¬ 
er  some  modifications  were  made  to  the  software;  these  are  discussed  in  Section 
II  of  this  report.  Detailed  updates  to  the  documentation  will  be  submitted  to  VELA 
Seismological  Center  under  separate  cover. 

Section  III  describes  work  performed  during  the  quarter  on  the  - 
•ALPA  evaluat.on  task.  Included  are  a  brief  description  of  routine  event  process  - 
ing,  results  of  a  study  to  separate  interfering  events  with  a  beamsteer  processor, 
a  amparison  of  the  signal  degradation  of  off-line  and  on-line  beamsteers,  and 
results  of  an  evaluation  of  two-component  Rayleigh  wave  processing. 

Section  IV  describes  work  performed  on  the  NORSAR  long  period 
evaluation  task.  The  evaluation  effort  during  the  quarter  was  directed  toward 


obtaining  and  processing  the  first  suite  of  events  and  noise  samples.  Due  to  pro¬ 
blems  encountered  in  finding  useable  data,  only  a  limited  number  of  events  and 
noise  samples  were  analyzed.  These  problems  are  discussed  along  with  the  gen¬ 
eral  data  quality.  Analysis  results  included  are  signal  analysis,  a  discussion  of 
the  ambient  noise  field,  some  initial  results  from  the  array  *  ’•ocessing  task,  and 
first  results  of  a  chirp  filter  study. 

Section  V  describes  work  performed  under  the  NORSAR  short  per¬ 
iod  evaluation  task.  Work  was  concentrated  on  analysis  of  the  general  quality  of 
the  132  channel  NDPC  data  tapes,  noise  analysis,  a  study  of  the  signal  character¬ 
istics,  and  an  investigation  of  array  beamforming  techniques. 

Section  VI  discusses  the  work  performed  on  the  LPE  task.  To 
date  the  LPE  evaluation  was  limited  by  the  quality  of  the  data  received  at  SAAC. 
However,  a  limited  amount  of  data  from  early  1971  have  been  processed  and  ana¬ 
lyzed.  Analysis  of  the  noise  field  for  this  time  period,  and  a  very  preliminary 
discussion  of  the  detection  capabiltiy  of  each  station  and  the  network  are  included. 
Also,  the  behavior  of  the  Mg-m^  classification  parameter  is  described. 

No  new  results  for  the  adaptive  processing  studies  are  available. 
Also,  no  results  from  the  spectral  analysis  task  are  available. 
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II. 


SOFTWARE  UPGRADES 


A.  Introduction 

During  the  past  quarter  a  few  modifications  were  made  to  the 
analysis  software.  Software  changes  were  made  to  EVPLOT  of  the  LP  array 
evaluation  package,  LXMERG  andLXDALL  of  the  LPE  package,  and  NSPED  of 
the  NORSAR  SP  package.  This  section  describes  the  program  modifications; 
appropriate  updates  to  the  documentation  will  be  submitted  to  VSC  under  separate 
cover  in  the  near  future. 

B.  Long-Period  Array  Evaluation  Package 
1.  EVPLOT 

During  the  past  quarter,  two  changes  were  made  in 
EVPLOT.  The  event  time  drawn  on  the  plot  is  now  updated 
by  the  time  interval  of  the  skipped  initial  segments,  if  any. 

An  amplitude  conversion  constant  and  alternative  rotation 
subroutine  call  are  provided  to  allow  the  program  to  be  used 
for  data  from  arrays  other  than  ALPA. 

C.  Long-Period  Experiment  Evaluation  Package 
1 .  LXMERG 

Several  changes  have  been  made  in  LXMERG  during  the 
past  quarter  on  the  basis  of  further  operating  experience. 

The  intermediate  tape  specification  card  now  uses  the  as¬ 
signed  tape  system  numbers,  rather  than  a  purely  internal 
tape  index  number.  A  console  interrupt  check  has  been  in  - 
eluded  in  the  main  merge  record  processing  loop  to  permit 
truncation  of  the  merge  interval  while  the  program  is  running. 
This  capability  is  valuable  for  partial  recovery  when  serious 
input  tape  problems  develop  late  in  the  merge. 


Input  tapes  have  been  found  to  contain  occasional  records 
with  spurious  timing  words.  To  permit  the  merge  program 
to  access  correct  data  subsequent  to  such  records,  an  add¬ 
itional  timing  test  has  been  included:  if  the  input  record 
time  is  greater  than  the  end  of  the  merge  interval,  the  record 
is  regarded  as  an  error  condition  and  is  ignored. 

Subroutine  TPROC  has  been  coded  and  is  now  used  for 
checking  the  timing  codes  in  the  input  intermediate  format 
merge  records.  The  complex  subscripting  involved  in  this 
checking  was  not  efficient  in  FORTRAN;  the  new  routine  is 
coded  in  ALC. 

Error  recovery  is  now  attempted  for  all  output  tape  error 
conditions.  Erase-gaps  and  record  repeats  are  attempted  up 
to  50  times.  On  the  final  library  tape,  a  maximum  of  11000 
output  records  are  written  on  each  tape,  since  the  system 
apparently  does  not  always  recognize  the  end-of-tape  reflec-  • 
tive  marker.  Finally,  the  start  and  stop  times  of  the  indivi¬ 
dual  output  tape  reels  are  printed  to  provide  more  convenient 
record  keeping. 

2.  LXDALL 

During  the  past  quarter,  minor  changes  were  made  in  the 
output  tape  format  to  conform  to  revised  specifications.  Out¬ 
put  error  recovery  coding  was  included:  such  conditions  as 
write-parity  now  cause  erase-gaps  and  output  re-attempts  to 
be  generated  up  to  fifty  times.  Changes  were  also  made  in 
the  tape  handling  conventions  to  allow  multiple  input  and  out¬ 
put  reels  to  be  processed  on  the  same  run. 
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D.  NORSAR  Short-Period  Array  Evaluation  Package 
1.  NSPED 

During  the  past  quarter,  initial  processing  of  large 
numbers  of  NDPC  tapes  was  carried  out,  and  several  pro¬ 
gram  changes  were  made  as  a  result.  It  was  assumed,  dur¬ 
ing  program  development,  that  the  NDPC  tapes  would  each 
contain  data  for  a  single  time  block  of  no  more  than  2,  5 
hours  duration.  In  practice,  each  tape  contains  data  for  two 
or  more  time  blocks  of  usussflly  30  minutes  each,  separated 
by  intervals  of  days,  and  not  always  in  correct  temporal  or¬ 
der.  Consequently,  the  tape  input  subroutine  READND  had 
to  be  modified  to  first  rewind  the  tape,  then  search  forward 
regardless  of  record  timing.  To  avoid  excessive  search  time, 
a  pause  and  interrupt  check  were  included  in  NSPED  to  allow 
the  analyst  to  instruct  the  program  at  the  beginning  of  each 
edit  search  whether  to  process  the  current  input  tape  or  re¬ 
quest  the  next  tape.  Changes  were  also  made  to  facilitate  in¬ 
put  error  recovery:  retry  attempts  are  now  made  up  to  thirty 
times  for  such  conditions  as  read-parity  and  short  records. 
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in.  ALPA  EVALUATION 

A.  Introduction 

Routine  processing  of  events  received  at  ALPA  has  continued 
during  this  quarter.  A  brief  description  of  this  work  appears  in  part  B.  In 
addition  several  investigative  efforts  have  been  performed.  Part  C  describes 
an  attempt  to  separate  interfering  events  with  a  beamsteer  processor.  In  part 
D  a  comparison  of  the  signal  degradation  of  off-line  and  on-line  beamsteers  is 
discussed.  Part  E  gives  recent  results  of  the  two-component  Rayleigh-wave 
processor  evaluation.  Future  plans  are  discussed  in  part  F. 

B.  Routine  Processing 

Approximately  fifty  events  which  occured  during  the  month 
of  June  have  been  processed.  This  includes  beamsteering,  bandpass  filtering, 
matched  filtering  when  necessary,  and  the  measurement  of  classification  para¬ 
meters.  Shear  wave  amplitudes  are  now  included  routinely  in  this  last  category. 
In  addition  to  the  above,  seven  events  from  known  underground  test  sites  have 
been  processed.  In  none  of  the  cases  considered  was  the  full  array  operational. 
Current  processing  is  concentrated  on  the  month  of  July.  It  is  believed  that 
data  from  the  latter  part  of  July  will  yield  a  higher  percentage  of  usable  sites. 

The  optimum  bandpass  filter  for  shear-wave  extraction  has 
been  determined.  A  suite  of  nine  events  was  selected  on  the  basis  of  having 
clear  arrivals  in  the  shear-wave  beamsteer.  Bodywave  magnitudes  ranged  from 
4.  7  to  6.  7,  and  deltas  ranged  from  27  degrees  to  77  degrees.  In  each  case  the 
component  having  the  largest  shear-wave  amplitude,  usually  the  transverse,  was 
bandpass  filtered  with  a  suite  of  filters  spanning  the  range  0.  012  Hz  to  0.  102  Hz. 
The  figure  of  merit  used  was  the  ratio  of  peak  signal  amplitude  to  RMS  noise 
value  in  the  bandpassed  beam.  In  five  of  the  nine  cases  the  best  performance 
was  achieved  with  an  0.027  Hz  to  0.  059  Hz  bandwidth;  in  three  cases  an  0.  020  Hz 
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to  0.051  Hz  bandwidth  was  best.  It  has  been  determined  previously  to  routinely 
use  an  0.025  Hz  to  0.055  Hz  bandpass  filter  on  the  surface-wave  beams.  In 
view  of  these  results  that  filter  also  appears  to  be  appropriate  for  the  shear- 
wave  beams.  The  degradation  in  signal  amplitude  resutling  from  bandpass 
filtering  was  computed  for  the  eight  cases  mentioned  above.  The  average 
degradation  was  found  to  be  1.7  dB. 

Seven  one-hour  noise  samples  were  edited  to  enable  com¬ 
parison  of  the  ALPA  and  NORSAR  noise  levels.  Figure  III- 1  shows  the  RMS 
values  of  these  samples  for  the  band  0.023  Hz  to  0.051  Hz,  averaged  across 
the  array.  Only  the  vertical  component  for  NORSAR  is  shown.  The  presence 
of  highly  site -dependent  very  long-period  noise  on  the  NORSAR  horizontal  seis¬ 
mometers  makes  an  array  average  of  questionable  value.  Based  on  this  brief 
sample  it  appears  that  the  vertical  component  noise  levels  at  the  two  arrays  are 
fairly  similar.  The  vertical  component  noise  at  NORSAR  in  this  band  is  largely 
dominated  by  the  micso seisms.  Thus  the  high  value  of  twelve  millimicrons  on 
day  191  is  due  to  the  presence  of  relatively  strong  microseismic  activity.  Note 
that  the  NORSAR  value  was  not  available  for  day  213. 

C.  Interfering  Event  Analysis 

On  26  July,  1971  the  SAAC  bulletin  reported  four  Solomon 
Islands  events  and  two  Sinkiang  events  occuring  between  0100  and  0300  GMT. 
The  Solomon  bodywave  magnitudes  given  were  6.  0,  5.  1,  6.  5,  and  6.  4;  the 
Sinkiang  magnitudes  were  6.  1  and  5. 1.  It  was  suggested  that  these  events  be 
processed  as  an  interfering  event  study.  In  an  effort  to  separate  the  events 
temporally  and  directionally  the  data  in  this  interval  were  broken  into  non- 

9 

overlapping  twelve  minute  and  forty-eight  second  time  gates.  High-resolution 
wavenumber  spectra  were  computed  for  each  time  gate  at  0.039  Hz  and  at  0.047 
Hz.  The  first  five  such  gates  which  spanned  the  surface  wave  arrivals  from  the 
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first  three  Solomon  events  and  the  first  Sinkiang  event  showed  peaks  only  in  the 
general  direction  of  the  Solomon  events.  The  sixth  gate  included  the  latter 
portion*,  of  the  fourth  Solomon  event,  the  second  Sinkiang  event,  and  an  mb=  6.  2 
Solomon  event  published  later  in  the  PDE  bulletin.  These  wavenumber  spectra 
showed  a  peak  at  3.  3  km/ sec,  345  degrees,  about  20°  from  the  correct  Sinkiang 
great  circle  azimuth.  This  peak,  while  quite  weak,  may  actually  represent  the 
Sinkiang  event. 

Surface  wave  beamsteers  were  formed  for  both  azimuths  us¬ 
ing  17  available  sites.  At  this  time  it  was  determined  that  the  first  Solomon 
Islands  event  has  a  surface  wave  magnitude  of  7.  7.  This  abnormally  high 
value  was  later  confirmed  by  the  PDE  bulletin  which  reported  Ms  =  7.  9.  As  a 
result  all  of  the  succeeding  Solomon  Islands  events  were  obscurred  by  the  coda 
of  this  event  and  could  not  be  detected.  The  Sinkiang  beam  significantly  atten¬ 
uated  this  event  but  not  enough  to  permit  detection  of  either  Sinkiang  event.  It 
is  likely  that  the  largest  of  the  Sinkiang  events  was  at  least  1.  5  surface-wave 
magnitude  units  smaller  than  the  obscurring  Solomon  Islands  event.  As  a  last 
resort  the  Sinkiang  beam  was  match  filtered  with  a  previously  recorded  Sinkiang 
master  waveform.  This  too  was  unsuccessufl  in  uncovering  either  Sinkiang 
event.  Apparently  the  differences  in  surface  wave  magnitudes  were  too 
great  to  permit  extraction  of  the  interfered  Sinkiang  events. 

D.  Comparison  of  On-Line  and  Off-Line  Beams 

The  purpose  of  this  study  was  to  compare  surface  wave 
beamsteer  outputs  formed  by  the  on-line  processor  with  those  formed  using 
off-line  software.  Differences  in  the  two  outputs  are  beamsteer  azimuths  and  site 
selction  prior  to  beam  formation.  Beamsteer  azimuths  used  in  the  off-line  soft¬ 
ware  correspond  to  the  expected  great  circle  path.  In  the  on-line  processor; 
however,  surface  wave  beamsteer  azimuths  are  preset  at  initialization  time. 
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Normally,  si;:  azimuths  in  increments  of  60°  have  been  chosen  for  on-line 
beamsteering.  Site  selection  for  off-line  beamsteers  are  made  by  the  analyst 
using  Calcomp  plots  of  each  site.  The  on-line  beams  are  formed  from  all  sites 
that  have  passed  the  software  quality  checks  for  a  given  15  second  block  of  data 
(Documentation  of  On-Line  Software,  1971).  It  is  therefore  possible  for  the  num¬ 
ber  of  sites  used  in  beamsteering  an  event  with  the  on-line  processor  to  vary  with¬ 
in  the  duration  of  the  event. 

Data  used'for  this  study  were  chosen  from  events  already  pro¬ 
cessed  through  the  off-line  software  and  reported  in  the  Final  Report  (Texas 
Instruments,  1971).  The  separation  between  great  circle  path  and  available  on¬ 
line  beamsteer  azimuths  for  these  events  ranged  from  zero  to  35  degrees.  T^e 
suite  included  both  high  and  low  signal -to -noise  ratio  events.  The  ten  events 
chosen  are  given  in  Table  III — 1 .  An  additional  five  events  were  selected  for  the 
comparison  but  were  not  accessible  by  the  on-line  plot  software. 

A  software  program  to  generate  Calcomp  plots  of  on-line  beams 
recorded  by  the  IIRSPS  system  was  borrowed  from  SAAC.  However,  on-line 
library  tapes  containing  the  desired  events  were  generated  by  the  Interim  ALPA 
system.  This  necessitated  considerable  software  change.  Calcomp  plots  of  the 
off-line  beamsteers  were  already  available  for  the  comparison. 

Signal  degradation  in  the  beamsteer  output  due  to  deviations  in  the 
beam  azimuth  can  be  anticipated  by  examinimg  the  beamsteer  response  in  f-k- 
space.  The  beam  pattern  for  the  seven  element  array  is  such  that  no  more  than 
3  dB  signal  degradation  would  be  expected  for  azimuthal  variations  as  large  as  40 
degrees.  Less  than  1  dB  degradation  is  espected  for  variations  less  than  20  de¬ 
grees.  Although  no  actual  comparison  of  on-line  and  off-line  beams  for  a  full 
19  element  array  is  possible  at  this  time,  it  should  be  noted  that  the  expected 
signal  degradation  for  off  azimuth  beams  is  more  severe.  In  order  to  be  certain 
that  no  more  than  1  dB  attenuation  in  signal  occurs,  an  azimuth  within  15°  of  the 
great  circle  path  is  required.  Likewise,  azimuthal  variations  of  35°  would 
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be  expected  to  cause  at  least  6  dB  signal  attenuation.  These  results  assume 

’  1  )  .  j 

that  the  event  propagates  across  the  array  as  a  plane  wave.  Practical  depart¬ 
ures  from  this  idealized  assumption  may  lead  to  results  inconsistent  with  this 

.  ,  1  .  '  1  ’  !  ! 

type  of  analysis.. 

i  '  , 

1  Table  III- 1  gives  the  great  circle  and  on-line  beam. azimuths 

\  } 

along  with  tlie  peak  amplitude  attenuation  measurements  for  each  of  the  events. 

1  '  l 

i 

The  numbers  given  are  the  radios  of  off-line  beam  peak  amplitude  to  on-line 
beam  peak  (amplitude,  expressed  in  dB.  In  all  cases  excepting  KUR/219/  13AL, 
the  differences  tietween  the  two  beam's  exceeded  two  dB  only  where  the  signal- 
to-noise  ratio  was  low.  In  other  words,  when  the  Signal  was  large  enough  to  be 

*  i 

clearly  measureable  the  diffei*ence  between  the  two  bearhs  was  two  dB  or  less. 

In  the  case  of  the  exception,  KUR/219/13AL,  the  on-line  beam  shows  a  large 
pul|Se  near  the  onset  of  the  vertical  component  Rayleigh -wave.  If  the  beam 

1  1  s  1 

measurement  is  made  on  this  pulse., the  on-lihe  beam  appears  4.  3  dB  larger  than 

1  i  i 

the  off-line  beam,.  It  is  possible,  however,  that  this  pulse  is  anomalous. 

'|i  .  ; 

1  '  I  i 

In  some  cases  when  the  separation  between  .great  circle  and 
on-line  beam  azimuths  is  large  some  crosstalk  between  the  horizontal  components 
is  observed,  i.  e.  the  Love-wave  appears  on  ,the  radial  component  or  vice  versa. 

If  the  LoVe-wave  and  horizontal  Rayleigh-wave  have  about  the  same  amplitudes 
this  does  not  appear  noticeable  in  the  amplitude  measurements,  but  might  be 
significant  in  the  case  of  a  large  Love-wave. 

1  i  '  1 

There  is  little  reason  to  believe  that  departure  of  the  on-line 
i  1 

beams  from, the  great  circle  azimuth  will  consistently  degrade  the  noise  suppress - 

1  1  '  ’ 

ion  of  the  beam.  In  light  of  the  results  above.it  would  appear  that  the  practice 
of  spacing  on-line  beams  at  60°  intervals  results  in  adequate  coverage  with  the 
subarray  considered.  This  is  particularly  true  of  ithe  vertical  component  where 
inter-mode  crosstalk  does  not  occur.  A  similar  comparison  with  the  full  array 
will  be  made  when  appropriate  data  become  available. 
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E.  Two-Component  Rayleigh  Wave  Processing 

Two -component  Rayleigh  wave  processing  consists  of  phase- 
shifting  the  in-line  horizontal  beam  by  90°,  applying  a  scale  factor  based  on  the 
ratio  of  vertical -to -horizontal  RMS  signal  levels,  and  stacking  the  vertical  and 
equalized  horizontal  beams.  In  this  study  6  events  and  19  noise  samples 
(Table  III- 2)  were  processed  over  4  bandwidths  and  two  scale  factors  to  deter¬ 
mine  the  noise  suppression  and  signal  enhancement  achieved  by  two-component 
processing.  The  bandwidths  used  were  the  following: 

0.020  Hz  to  0.060  Hz 
0.020  Hz  to  0.047  Hz 
0.027  Hz  to  0.060  Hz 
0.027  Hz  to  0.047  Hz 

The  two  horizontal  scale  factors  were  1.0  and  1.2. 


The  measure  of  the  average  increased  noise  suppression  was 
obtained  by  comparing  the  integrated  power  spectra  (at  each  bandwidth  and  scale 
factor)  of  the  two-component  processed  noise  with  that  of  the  beamsteered  vert¬ 
ical  component.  The  power  values  of  these  differences  in  dB  are  given  in  Tabic 
III- 3. 


In  the  case  where  the  vertical  and  in-line  horizontal  beams  had 
the  same  noise  power  density  and  were  incoherent,  noise  rejection  would  be  1.  3 
dB  (Vn~  ).  On  the  average,  two-component  processing  with  a  scale  factor  of 
1. 2  yields  only  slightly  greater  noise  rejection  than  the  1. 3  dB  figure,  while  a 
scale  factor  of  1.0  gives  from  0.7  to  1.0  dB  greater  rejection  than  the  1.  3  dB 
figure.  Perhaps  even  more  significant  is  the  average  noise  suppression  obtained 
by  two-component  processing  relative  to  the  smellier  of  the  vertical  and  in-line 
Rayleigh  wave  components  (Table  IH-4).  The  average  noise  rejection  using  a 
scale  factor  of  1.0  is  almost  one  dB  greater  than  using  a  scale  factor  of  1.2,  at 
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Associated 

Event 

Name 

No. 

of 

Sites 

Date/ 
Hour  of 
Event 

Date/ 

Hour  of 
Noise 
Sample 

Ms 

mb 

IRA/ 242/ 1&*  1 

. 6  “ 

242/1600 

242/1500 

4.  1 

5. 1 

CHI/220/11* 

5 

220/1100 

220/1000 

4.0 

4.7 

IRA/211/13 

6 

211/1300 

211/1000 

— 

H 

CRS/287/06* 

7 

287/0600 

287/1300 

4.6 

KYU/207/07* 

8 

207/0700 

207/0600 

5.8 

6. 1 

GRC/184/00 

6 

184/0000 

184/0300 

3.9 

5. 1 

CHI/212/ 1 3* 

6 

212/1300 

212/1200 

4.  5 

5.  5 

KAM/242/00* 

7 

242/0000 

241/2300 

4.8 

5.2 

KAM/231/Q9 

9 

231/0900 

234/1000 

3.  1 

4.7 

KAM/ 285/ 22 

7 

225/2200 

225/2100 

3.  3 

4.6 

BUR/225/07 

9 

225/0700 

225/0300 

3.7 

4.7 

KAM/242/03 

6 

242/0300 

241/2200 

3.6 

4.8 

KAM/242/07 

9 

242/0700 

242/0600 

2.6 

3.8 

CAU/243/04 

7 

243/0400 

243/1200 

3.0 

ALB/241/10 

9 

241/1000 

241/08CO 

3.  3 

4.4 

SIN/ 240/05 

6 

240/0500 

240/0000 

— 

— 

WKZ/357/07 

7 

357/0700 

357/0600 

4.  3 

6.  1 

EKZ/351/07 

6 

351/0700 

351/0600 

3.  3 

5.  3 

EKZ/147/07 

3 

147/0700 

147/0300 

*— 

— 

*  Events  used  in  the  Bignal  enhancement  study 


TABLE  III-2 

19  EVENTS  AND  ASSOCIATED  NOISE  SAMPLES  USED  IN  EVALUATION  OF 

TWO-COMPONENT  PROCESSING 
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Improvement  Of  T wo -Component  Over  Vertical  Trace  (dB) 

0.020  to  0.  020  to  0.  027  to  0.  027  to 

Associated  0.060  Hz  0.047  Hz  0.060  Hz  0.047  Hz 

Event  Scale  Factor  Scale  Factor  Scale  Factor  Scale  Fact 
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NOISE  SUPPRESSION  ACHIEVED  BY  TWO-COMPONENT  PROCESSING 


Bandwidth 

(Hz) 

Scale  factor 
for  the 
Horizontal 
Component 

Average 

Woise 

Suppression 

(dB) 

0.020-0.060 

1.0 

1.9 

0.020-0.060 

1.  2 

1.0 

0.020-0.047 

1.0 

1.6 

0.020-0.047 

1.  2 

0.8 

0.027-0.060 

1.0 

2.0 

0.027-0.060 

1.2 

1. 1 

0.027-0.047 

1.0 

1.6 

0.027-0.047 

1.  2 

0.7 

Difference 
In  Achieved 
Noise  Suppression 
(Scale  factors  1.  2) 
-(Scale  factor=  1.0) 


TABLE  III -4 

AVERAGE  TWO-COMPONENT  NOISE  REJECTION  IMPROVEMENT  RELATIVE 
TO  THE  SMALLER  OF  VERTICAL  AND  HORIZONTAL  COMPONENTS 


all  bandpasses.  The  scale  factor  of  1.0  achieves  noise  suppression  of  from  0.  3 
to  0.  7  dB  beyond  the  level  expected  for  incoherent  noise  where  the  input  noise 
level  equals  the  lower  of  the  vertical  and  in-line  horizontal  noise  levels. 

The  signal  gains  (Table  III- 5)  were  obtained  by  taking  the  ratio  of 
the  maximum  peak  amplitude  of  the  wideband  processed  two -component  output 
and  the  amplitude  of  the  wideband  vertical  Rayleigh  wave  at  the  same  frequency. 
The  signal-to-noise  ratio  improvement  (Table  III-5)  was  calculated  as  the  differ¬ 
ence  in  peak  signal  amplitude  over  RMS  noise  between  the  two-component  pro¬ 
cessor  output  and  the  vertical  Rayleigh  wave.  The  two-component  processor  is 
slightly  more  effective  when  the  bandpass  is  extended  to  0.06  Hz  on  the  high 
frequency  end  than  when  it  is  cut  off  at  0.047  Hz.  This  improved  signal-to-noise 
ratio  enhancement  results  from  greater  noise  suppression  at  the  higher  fre¬ 
quencies  as  seen  by  the  averages  of  Table  HI-3.  In  most  cases  the  microseismic 
noise  at  ALPA  propagates  as  surface  waves  coming  from  the  South.  Since  the 
signals  of  Table  HI-5  tend  to  come  from  the  North,  the  processor  can  reject 
the  microseismic  Rayleigh-wave  noise  by  more  than  the  1.  3  dB(V?Tj  figure. 

This  capability  is  noted  when  the  bandpass  is  extended  to  include  the  micro - 
seisms. 

Scaling  the  horizontal  trace  by  1.0  yields  slightly  better  results 
than  the  use  of  the  1.  2  scalar.  Apparently  increasing  the  value  to  1.  2  increases 
the  noise  level  more  than  the  signal  level  in  the  processor  output. 

These  results  indicate  that  the  0.025  to  0.055  Hi  bandpass  pre¬ 
viously  selected  for  array  processing  and  matched  filtering  is  also  suitable  for 
two-component  processing.  The  two-component  processor  achieves  about  two 

dB  signal -to -noise  ratio  enhancement  over  the  bandpassed  vertical  component 
beam. 
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SIGNAL  GAINS  AND  SIGNAL -TO -NOISE  RATIO  ENHANCEMENT 
FOR  TWO -COMPONENT  RAYLEIGH  WAVE  PROCESSING 


F.  FUTURE  PLANS 

Routine  processing  of  events  with  a  two  to  three  month  lag  between 
the  time  of  occurence  and  the  processing  of  the  events.  This  lag  ensures  that 
PDE  information  is  available  prior  to  processing  of  the  data.  When  the  full 
array  data  become  available  the  events  will  be  processed  with  all  available  sites 
and  with  a  seven  site  subarray.  This  will  permit  assessment  of  the  relative 
capabilities  of  the  two  arrays.  As  events  are  processed,  those  suitable  for  use 
as  master  waveforms  will  be  retained  for  this  purpose. 

Routine  analysis  of  noise  samples  spaced  at  ten  day  intervals  will 
continue.  It  is  contemplated  that  this  can  be  done  most  efficiently  in  batches 
covering  two  or  three  month  intervals. 
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IV. 


NORSAR  LONG-PERIOD  EVALUATION 


A.  Introduction 

The  evaluation  effort  during  the  second  quarter  has  been  directed 
toward  obtaining  and  processing  the  first  large  suite  of  events  and  noise  samples. 
Fourteen  events  were  edited  and  are  now  in  analysis;  event  proce,  jing  almost 
has  been  completed.  Also,  eleven  noise  samples,  including  two  long-duration 
(  >3  hours)  samples,  have  been  obtained  and  processing  is  almost  complete. 

A  larger  set  of  events  would  have  been  processed  except  that  an  unusually  large 
amount  of  time  was  spent  finding  and  editing  suitable  data.  The  NORSAR  array 
was  just  becoming  operational  during  this  period. 

Results  of  the  evaluation  to  date  are  presented  in  the  following 
sections.  Section  B  reviews  data  quality;  section  C  discusses  signal  analysis; 
section  D  discusses  the  ambient  noise  field;  section  E  discusses  some  initial  re¬ 
sults  from  the  array  processing  task;  section  F  presents  the  first  results  of 
chirp  filtering  the  array  output  beams.  No  information  on  NORSAR  long-period 
detection  capability  or  discrimination  capability  is  yet  available;  preliminary  data 
on  these  topics  will  be  obtained  in  the  near  future. 

B.  General  Data  Quality 

The  small  amount  of  data  processed  during  this  quarter  reflects 
the  difficulty  encountered  in  obtaining  good  data.  This  difficulty  is  primarily  the 
result  of  non-continuity  in  time  of  the  data  (mostly  ^caused  by  transmission  pro¬ 
blems)  CtlAvI  incomplete  event  information. 

Events  edited  during  this  quarter  were  obtained  from  the  interval 
day  168  to  day  216  (17  June  to  4  August,  1971).  Since  available  "PDE"  bulletins 
only  covered  the  first  part  of  this  period  the  SAAC  bulletins  were  the  primary 
source  of  event  information.  The  total  number  of  events  from  the  Eurasian 
landmass  (including  Russia  but  excluding  East  and  West  Europe)  given  by  the 
SAAC  bulletin  was  123. 
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Since  this  was  preliminary  analysis,  we  chose  only  those  events 
which  were  either  visible  on  the  raw  data  or  which  were  from  areas  of  interest 
and  listed  by  PDE.  Of  these  123  events  only  nine  contained  both  visible  surface 
waves  and  reasonably  continuous  time  data.  Actually,  twelve  events  were  pro¬ 
cessed  since  three  others  from  areas  of  interest  were  listed  by  PDE  .  A 
breakdown  of  the  reasons  these  114  events  were  initially  rejected  is  as  follows: 

•  data  was  not  available  for  43  events  (NDPC  data  was  not  present 
on  the  SAAC  library  tapes  for  the  period  from  day  174  to  day  183) 

•  28  events  had  large  time  gaps  in  the  data  which  precluded  ana¬ 
lysis 

•  4  events  were  definitely  interfered  with  by  other  events  (interfering 
events  listed  by  SAAC) 

•  14  events  appeared  to  be  interfered  with  by  other  events 

•  5  events  had  uncorrectable  amplitude  spikes  occuring  in  the  sur¬ 
face  wavetrain 

•  13  events  were  not  visible  (three  were  subsequently  confirmed  by 
PDE  and  were  then  used) 

•  5  events  were  mislocated  by  SAAC 

•  1  event  was  too  deep  to  be  of  interest 

•  1  event  occurred  when  the  data  on  all  sites  and  components  was 
identical 

Almost  all  of  the  data  not  rejected  due  to  time  gaps  contained  some 
of  these  gaps.  When  the  gaps  were  not  excessive  or  if  they  were  not  present  in 
the  surface  wave  trains,  the  data  were  processed  with  allowances  made  for  the 
missing  time. 
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Contrary  to  our  original  suspicions,  data  dropout  now  appears  to 
be  mostly  random  and  unrelated  to  any  seismic  wave  arrivals  at  NORSAR.  In¬ 
formation  on  the  dropout  rate  is  now  available  but  has  not  yet  been  tabulated. 

The  ten  day  loss  of  data  apparently  was  due  to  an  undetected  transmission  failure. 
It  is  believed  that  operational  procedures  at  SAAC  have  been  revised  to  prevent 
a  recurrence  of  this  type  of  data  loss.  Another  source  of  data  loss  is  use  of 
the  TAL  for  vocal  communication  with  NDPC  (data  cannot  be  transmitted  when 
the  modem  is  in  voice  mode).  The  amount  of  data  lost  from  this  cause  has  not 
been  investigated;  however,  communication  to  NDPC  via  the  computer  console 
typewriter  should  be  used  whenever  possible  to  prevent  interruption  of  the  data. 

C.  Signal  Analysis 

During  this  quarter,  four  additional  events  with  large  S/N  ratio 
were  edited  for  signal  analysis.  The  six  events  presently  being  analyzed  are 
listed  in  Table  IV- 1. 

This  analysis  consists  of  determining  signal  characteristics  that 
will  affect  array  processing.  The  questions  to  be  answered  are: 

•  Are  the  signals  similar  at  each  site? 

•  Do  events  arrive  along  the  great  circle  route  and  at  the  expected 
times  ? 

•  What  is  the  spectral  content  of  the  signals  and  the  variation  by 
region? 

•  What  is  the  signal  duration? 

•  Are  there  multipathed  arrivals  and  what  are  the  secondary  arrival 
directions? 

Signal  similarities  have  been  measured  for  all  six  events.  Table 
IV- 2  presets  the  average  similarity  (correlation  coefficient)  Detween  the  reference 
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LIST  OF  EVENTS  FOR  SIGNAL  ANALYSIS 
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TABLE  IV- 2 


SIGNAL  SIMILARITY 


Event 

S  vertical 

S  radial 

LQ  trans¬ 
verse 

LR  vert¬ 
ical 

LR  radial 

TUR/161/09NL 

0.  840 

0.906 

0.916 

0.893 

0.  836 

TUR/143/01NL 

0.921 

0.  900 

0.867 

0.863 

0.810 

TUR/126/04NL 

0.  969 

0.923 

0.950 

0.885 

0.  836 

KUR/190/16NL 

0.813 

0.815 

0.  792 

KUR/213/02NL 

0.917 

0.930 

0.890 

0.884 

0.  869 

CHI/208/13NL 

0.907 

0.929 

0.  898 
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site  and  the  other  sites  for  the  vertical  and  radial  S,  transverse  Love  and  vert¬ 
ical  and  radial  Rayleigh  modes.  As  can  be  seen,  the  coefficients  are  above 

0.  80  indicating  good  similarity  across  the  array.  ’ 

» 

Power  spectra  were  also  generated  for  each  event,'  and  Figures 
IV-1,  IV-2  and  IV-3  show  spectra  from  three  different  areas.  The  locations  of 
the  three  Turkish  events  are  very  close.  TUR/126/04NL  and  TUR/161/09NL  are1 

i  i 

within  0.  1  degrees  of  each  other  and  their  spectra  are  nearly  identical.  The  ! 
spectra  for  TUR/143/01NL  (Figure  IV-1)  which  is  approximately  3.  5  degrees 
from  the  other  two,  is  essentially  the  same  as  the  others  although  the  width  of 
the  -6  dB  band  (6  dB  down  from  the  peak  power)  for  the  vertical  component  is 
narrower  (0.025  Hz  -  0.06  Hz  for  TUR/143/01NL  versus  0.025  Hz  -  0.070  Hz^ 
for  TUR/126/04NL). 

1  *  !  ' 

The  Kurile  events  are  separated  by  fourteen  degrees  however,  they 
both  peak  at  0.  040  Hz  on  the  verticals  and  the  widths  of  the  -6  dB  band  for  the 
verticals  are  the  same  (0.025  Hz  to  0.060  Hz). 

? 

The  Chinese  event  CHI/208/13NL  has  more  high  frequency  enerjgy 
than  the  other  events  (width  of  the  -6  dB  band  is  0.040  Hz  to  0.079  Hz)  and  the 
vertical  peak  is  0.  055  Hz. 

In  order  to  determine  the  optimum  passband  for  more  routine 

J  ' 

processing,  the  bandwidth  behavior  of  many  events  from  different  regions  must 
be  observed.  If  significant  signal  energy  is  present  at  higher  frequencies  !for 
events  from  some  regions,  passbands  tailored  for  each  region  may  be  necessary. 
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Power  (<J B  relative  to  1  at  0,04  Hz) 


Frequency  (Hz) 


FIGURE  IV- 2 

3  COMPONENT  POWER  SPECTRA  FOR  EVENT  KUR/190/16NL 

SITE  3 
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Power  i?®  relative  to  1  mjf2/H*  at  0.  04  Ha) 


FIGURE  i'V-3 

3  COMPONENT  POWER  SPECTRA  FOR  EVENT  CHI/208/  13NL 

SITE  3 
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D.  Noise  Analysis 

A  preliminary  analysis  of  the  NORSAR  long-period  noise  was 
conducted  using  ten  noise  samples  taken  at  approximate  ten  day  Intervals  during 
the  period  from  day  121  to  day  220  (May  1(  1971  to  August  8,  1971)  Tabic  IV-3 
is  a  list  of  the  noise  samples  and  their  lengths. 

A  typical  3-component  noise  power  density  spectrum  Is  shown  in 
Figure  IV-4.  The  system  response  is  not  removed,  but  the  spectrum  is  nor¬ 
malized  to  1  m^j2/Hz  at  0.04  Hz.  The  16-18  second  microseismic  peak  dominates 
the  spectrum  with  levels  at  30  to  40  dB.  There  is  generally  a  small  notch  be¬ 
tween  .  090  and  0.  10  Hz  (10-1 1  seconds).  The  three  components  are  generally 
within  3  dB  from  0.047  to  higher  frequencies. 

Below  the  microseismic  peak  long-period  noise  often  appears  on 
the  horizonthls  (Figure  IV -5).  The  long-period  horizontal  noise  is  highly  variable 
from  site  to  site  and  is  observed  in  the  signal  processing  band  (20-40  second 
periods).  Figure  IV-6  shows  the  variations  of  the  RMS  ground  motion  in  the 
20-40  second  period  band  as  a  function  of  site.  As  can  be  seen,  the  horizontal 
RMS  values  fluctuate  considerably  more  than  the  vertical  values  (9  dB  vs  3  dB). 

At  the  mlcroscsimic  peak  noise  levels  are  reasonably  similar  across  the  arrav* 
and  among  components  (Figure  IV-7).  To  date  no  frequency  reeponse  calibrations 
have  been  available,  however,  the  small  variation  in  microseismic  peak  noise 
level  indicates  reasonably  good  equalization  between  sites. 

The  20-40  second  array  average  RMS  noise  levels  were  plotted 
versus  time  to  investigate  time  variability  of  the  NORSAR  long-period  noise 
field,  (Figure  IV-8).  The  fluctuations  in  the  horizontal  noise  level  is  seen  to 
be  much  larger  than  the  vertical  noise  (12  dB  maximum),  but  no  long  term  trend 
is  indicated.  RMS  levels  in  this  band  ranged  from  5  to  12  mjj  for  the  vertical, 

6  to  20  mj|  for  the  North-South,  and  5  to  15  m^  for  the  East- West.  These 
values  are  similar  to  those  observed  at  ALPA. 
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TABLE  IV -3 
NORSAR  NOISE  SAMPLES 


DAY 

LENGTH 

121 

39m-24  sec 

135 

lhr-33m-52sec 

Ml 

lhr-l6m-48sec 

150 

57m-28scc 

164 

4hr»37m-20sec 

171 

lhr-16m-48scc 

191 

51m-12sec 

201 

59m-44sec 

21 1 

lhr-12m-32sec 

220 

4hr-24m-32sec 
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FIGURE  IV -5 

EXAMPLE  OF  STRONG  LONG  PERIOD  ENERGY  ON  HORIZONTALS 
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Figures  IV-9  and  IV-10  show  the  multiple  coherence  of  the  two 
four  hour  samples  at  day  164  and  day  220  for  all  three  components.  The  dots 
above  the  curves  indicate  the  peaks  in  the  power  density  spectrum  at  these  days. 
Figure  IV-9, for  day  164,  shows  the  two  noise  power  density  spectrum  peaks 
(at  16-18  seconds  and  8.  5  seconds)  to  be  fairly  coherent  on  the  vertical  (0.9  and 
0.  8  respectively)  but  only  moderately  coherent  on  the  horizontals  (  in  fact  the  East- 
West  component  does  not  have  an  8  second  period  peak).  The  noise  sample  shows 
low  coherence  at  periods  below  the  micro  seismic  peak  between  .07  to  0.  10  Hz. 

Figure  IV-10,  for  day  220,  shows  a  generally  broader  coherent 
peak  of  about  0.  9  on  all  three  components  between  0.  05  and  0.  08  Hz.  The  null 
at  0.07  to  0.  10  Hz  seen  on  day  164  is  not  present  here  (there  was  no  null  in  the 
power  density  spectra  for  this  day).  The  coherence  above  0.08  Hz  generally 
drops  with  increasing  frequency,  although  there  is  a  "plateau"  in  the  vertical 
coherence  between  0.  08  and  0.  15  Hz. 

These  coherence  results  are  not  encouraging  with  regard  to  multi¬ 
channel  filtering;  based  on  this  very  preliminary  analysis  it  does  not  appear  that 
multichannel  filtering  will  significantly  outperform  beamsteering  in  noise  rejec¬ 
tion. 

High  resolution  frequency-w’avenumber  spectra  were  calculated  for 
all  noise  samples  to  determine  direction  and  velocity  of  the  coherent  noise  power 
spectra  peaks.  Table  IV-4  shows  the  results  of  the  vertical  component  for  each 
noise  sample.  The  most  common  propagation  direction  was  90°  to  100°,  although 
noise  appeared  from  every  quadrant.  Since  the  signals  of  interest  will  have 
azimuths  from  0°  to  120°,  the  noise  and  signal  will  have  similar  azimuths  in  many 
cases  and  less  than  Vn~  array  processing  noise  rejection  might  be  expected  for 
these  cases. 
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MULTIPLE  COHERENCE  DAY  164  SITE  1  PREDICTED  FROM  SITES 
2,  3,  4,  7,  8,  9,11,13,  14,15,16, 

17, 18,  19,  22 
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TABLE  IV-4 

FREQUENCY,  VELOCITY  AND  AZIMUTH  OF  NOISE 

SAMPLES 


Day 


Frequency 


121 

135 

141 

150 

164 

171 

191 

201 

211 


.059 
.059 
.055 
.  066 
.059 
.051 
.051 
.059 
.059 
.059 


Velocity 

3.  6 
3.  8 
3.  8 

3.6 
3.  8 

3.7 
3.7 
3.7 
3.  8 


Azimuth 

197° 
92° 
108° 
5° 
240b 
195° 
90° 
60° 
100° 
30  3? 


220 


3.7 


i 


i 


! 


I 

»  • 

E.  Arra^  Processing  Pejrfbrmance 

i 

i  s 

i  This  analysis  is  directed  toward  comparing  the  performance  of 
donventioinal  beamsteering  and  multichannel  filtering  (MCF)  at  NORSAR.  The 
response  of  both  processors  to  signals  can  be  made  to  be  essentially  identical 

t  *  * 

and  so  their  perforinance  can  be  compared  in  terms  of  their  ability  to  suppress 

!  j  ' 

noise  in  the  output  beam. 

i  1  . 

1  '  Two  long  duration  and  one  medium  duration  noise  samples  have 

been  used  for  MCF  design  and  testing.  First,  a  portion  of  each  sample  was  used 
to  estimate  a  set  of  crpssppw^r  matrices,  then  p  vertical  component  MCF  was 

I  : 

designed  using  a  signal -to -noise  ratio  of  4.0  and  a  dispersive  signal  model  steer¬ 
ed  to  ian  azimuth  of  76°.  The  beamsteer  filter  and  the  MCF  were  applied  to  two 
sections  of  the  data.  The  first  section,  within  the  MCF  design  gate,  produced 
the  "on-design"  output,  and  the  second  section,1  not  used  in  the  design  gate,  pro¬ 
duced  the  "off-design"  output.  The  amount  of  data  foir  each  sample  used  for  this 

I  I  ! 

analysis  are  given  in  Table  IV-5. 

I 

!  Input  power  was  computed  by  averaging  the  vertical  component 

power  from  each  site  used  in  the  MCF  design.  The  MCF  and  beamsteer  filters 

I  !  ! 

utilized  the  same  sites.  The  ratios  of  the  average  input  power  to  the  output 
power  expressed  in  dB,  for  the  MCF  and  beamsteer  filter,  in  the  frequency  band 
of  0.023  to  0,051  Hz,  areigiVen  in  Table  IV-6.  These  results  are  similar  to  those 
obtained  from  ALPA.  Operating  in  the  desigh  gate,  the  MCF  achieves  5.  7,  2.  7, 

I 

and  0.  8  dB  of  additional  noise  suppression  over  the  beamsteer.  Outside  of  the 
design  gate,  the  relative  performance  drops  2.  3,  0.  5,  and  -1.  3  dB.  Note  that  in 
only  one  case  (the  MCF  processor  foir  day  164)  was  the  in  noise  rejection 

exceeded  for  off -design  noise. 

1  i  , 

'ii  _  i  i 

The  degradation  of  MCF  performance  for  off-design  noise  can  be 

due  to  either  insufficient  statistic^  or  non-stationarity  of  the  noise.  In  two  of  the 

i  i  ' 

1  1  1  i 
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TABLE  IV-5 

MCF  AND  BS  FILTER  DESIGN  PARAMETERS 


Sample 

Number 

Of  Sites 

Design 

Gate 

Duration 

(sec) 

On- 

Design 

Duration 

(sec) 

Ofi- 

Design 

Duration 

(sec) 

NOIS/164 

13 

12800 

5120 

5120 

16 

5632 

3072 

3072 

1 

15 

11520 

5888 

2816 

TABLE  IV-6 

ARRAY  PROCESSOR  NOISE  SUPPRESSION  * 


Noise  In  /  Noise  Out  (dB)  vertical  component 

Vn“ 

On-Design 

Off-Design 

Sample 

Gain 

BS 

MCF 

BS 

MCF 

NOIS/164 

11.1 

7.2 

12.9 

9.3 

11.6 

EKZ/145 

12.0 

7.6 

10.3 

5.8 

6.  3 

NOIS/220 

11.8 

10.0 

10.8 

9.2 

7.9 
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three  samples  over  three  hours  of  noise  were  used  to  estimate  the  spectral  mat¬ 
rices  for  considerably  fewer  than  the  full  22  sites  potentially  available  and  so 
the  degradation  in  off-line  performance  is  most  likely  due  to  the  non-stationarity 
of  the  noise  field.  Based  on  these  preliminary  results,  it  does  not  seem  likely 
that  MCF  processing  will  offer  significant  improvement  relative  to  beamsteering 
at  NORSAR.  This  task  will  be  continued,  however,  with  several  more  long 
noise  samples  so  that  a  more  definitive  answer  as  to  MCF  effectiveness  can  be 
obtained. 

F.  Matched  Filter  Performance 

Chirp  filter  performance  at  NORSAR  fens  been  measured  for  nine 
Eurasian  events  listed  in  Table  IV-7.  The  perform  assef  was  measured  by  taking 
the  ratio  of  the  maximum  amplitude  of  the  chirp  HH'«r  to  the  maximum  amplitude 
of  the  array  beam,  after  application  of  a  bandpass  filter  with  the  same  passband 
as  the  chirp.  This  measurement  is  valid  because  both  the  chirp  and  bandpass 
filter  will  have  the  same  RMS  noise  output  level.  At  the  present  time,  the  best 
passband  for  NORSAR  has  not  been  determined.  It  is  likely  that  the  band  will  be 
close  to  that  for  the  ALPA;  thus  a  .  025  to  .055  H*  band  was  selected  for  these 
events.  The  chirp  filter  length  was  chosen  for  each  event  after  rough  measure¬ 
ments  of  duration  and  dispersion  were  made  on  the  waveforms.  No  master  wave¬ 
form  filtering  has  been  done  yet;  however,  two  events  which  have  been  chosen  as 
masters  have  been  matched  against  themselves. 

The  chirp  filtering  results  to  date  are  shown  in  Table  IV-7  and  the 
master  events  matched  (and  bandpassed)  to  themselves  are  shown  in  Table  IV-8. 
The  following  observations  may  be  made  for  this  small  ensemble: 

•  The  Love  wave  gains  arc  generally  smaller  than  Rayleigh  wave 

gains  except  for  two  Kurile  events.  These  two  appear  to  have  the 
same  source  mechanism  since  their  waveforms  are  almost  identical 
point  for  point. 
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TABLE  IV -7 

CHIRP  FILTER  IMPROVEMENTS 


Improvement  Chirp/ BP  (dB) 

0.025  -  0.055  Hb 

Event 

mb 

LQT 

LRV 

LRR 

Mb 

KUR/190/16 

4.9 

69.9 

6.4 

3.9 

5.3 

4.8 

KUR/ 191/03 

4.8 

69.8 

6.1 

4.0 

3.9 

4.4 

KUR/191/09 

4.6 

69.2 

ND* 

ND 

ND 

— 

KUR/ 21 3/02 

5.7 

65.1 

3.6 

6.3 

5.6 

4.9 

SIB/212/01 

3.5 

55.3 

— 

2.5 

0.7 

3.4 

CHI/  208/ 13 

5.3 

45.8 

0.4 

1.3 

1.2 

4.6 

BSA/210/19 

4.5 

22.2 

— 

0.3 

— 

2.8 

5RS/200/20 

3.8 

19.8 

-1.6 

2.1 

3.3 

3.1 

WRS/191/17 

5.3 

20.  3 

ND 

ND 

ND 

— 

*  Nc<t  Detected 


TABLE  IV-8 

MASTER  APPLIED  TO  MASTER  IMPROVEMENTS 


Improvement  (dB)  0. 

025  -  0.055  Hb 

Event 

LQT 

LRV 

LRR 

KUR/190/16 

11.5 

6.6 

9.7 

CHI/ 208/ 13 

4.7 

3.5 

4.5 
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•  The  waveform  of  ♦he  China  event  was  not  well  dispersed  and  was 
of  unusually  she  rt  duration.  Chirping  would  not  be  expected  to 
be  effective  for  that  event, 

•  The  radial  and  vertical  component  Rayleigh  chirp  improvements 
are  close,  as  would  be  expected. 

•  The  chirp  filter  lengths  chosen  may  not  be  optimum  in  a  few 
cases  and  some  additional  improvement  may  be  available. 

•  The  two  e\  »nts  matched  against  themselves  gave  about  4  dB  more 
improvement  than  that  obtained  by  chirp  filtering. 

As  the  signal  and  noise  structure  become  better  defined  it  is  ex¬ 
pected  that  additional  chirp  improvement  will  be  achieved.  Some  of  the  above 
events  will  be  reprocessed  to  obtain  any  additional  gain  available. 

G.  Future  Plans 

During  the  third  quarter,  evaluation  will  be  directed  at  further 
definition  of  array  behavior  with  events  of  moderate  magnitude.  Definition  of 
passband  limits  for  event  processing  will  have  high  priority.  Reducing  the  num¬ 
ber  of  trial  bandpass  filters  used  reduces  proportionally  the  amount  of  processing 
time  needed. 

Noise  analysis  will  continue  with  routine  noise  sample  edits  taken 
at  approximately  ten  day  intervals  as  well  as  edits  of  at  least  one  long  duration 
sample  each  month  if  data  are  available. 

Chirp  filteri.  will  continue  and  matched  filtering  will  be  begun. 

A  certain  amount  of  reworking  of  the  present  suite  of  events  is  planned  through 
further  optimizing  of  the  chirp  results.  Matched  filtering  with  master  waveforms 
will  begin  and  standard  discriminants  will  be  computed. 
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The  present  goal  for  the  number  of  new  events  to  be  editted  and 
processed  for  the  third  quarter  is  at  least  45.  Whenever  possible,  the  PDE 
bulletin  will  be  used  for  event  information,  however  the  SAAC  bulletin  will  be 
used  for  the  lower  magnitude  events.  Events  with  magnitudes  between  4.  0  and 
5.  5  will  be  of  primary  interest. 
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V. 


NORSAR  SHORT -PERIOD  EVALUATION 


The  NORSAR  short-period  array  is  a  100 -kilometer -diameter 
array  with  22  subarrays  each  with  6  sensors.  Every  site  is  located  within  a 
different  subarray  of  the  short-period  array.  Each  short -period  subarray  con¬ 
tains  a  central  sensor  surrounded  by  a  5- sensor  ring  appriximately  7  kilometers 
in  diameter.  In  this  section  subarray  1  refers  to  subarray  1A,  subarrays  2 
through  8  to  subarrays  IB  through  7B,  and  subarrays  9  through  22  to  subarrays 
1C  through  14C.  Within  a  subarray,  sensors  0  through  4  refer  to  sensors  in 
the  surrounding  ring,  starting  at  approximately  0°  relative  to  the  central  sensor 
of  the  subarray.  Sensor  5  refers  to  the  central  sensor  of  the  subarray.  Thus 
sensor  0  of  subarray  1  is  sensor  01A01  in  the  official  nomenclautre,  and  sensor 
5  of  subarray  22  is  sensor  14C00. 

NORSAR  short-period  array  analysis  effort  to  date  has  been  con¬ 
centrated  in  four  areas: 

•  General  quality  of  the  132  channel  data 

•  Noise  analysis 

•  Study  of  signal  characteristics 

•  Determination  of  array  beamforming  techniques 

Subsection  A  briefly  covers  the  quality  of  data  transcription  for  the 
NORWAY  NDPC  tapes.  Subsection  B  describes  highlights  of  noise  analysis  con¬ 
ducted  on  the  14  noise  samples  processed.  Subsection  C  presents  the  results  of 
signal  analysis  performed  on  selected  events.  Subsection  D  discusses  attempts 
to  improve  beamforming  performance  through  the  determination  of  a  set  of  adjusted 
propagation  delays  for  different  event  epicenters  and  through  the  selective  weight¬ 
ing  of  subarrays  during  array  summation  (using  a  diversity- stack  technique). 


V-l 


A. 


Data  Transcription  Quality 


Thus  far,  the  NORSAR  short-period  edit  program  has  attempted  to 
transcribe  22  signals  and  a  5-minute  noise  sample  preceding  each  signal  from 
the  NORWAY  NDPC  tapes.  Fifteen  signals  and  14  noise  samples  were  success¬ 
fully  edited.  Table  V-l  shows  epicenter  data  for  the  15  signals.  No  noise  sample 
was  obtained  for  Event  3.  About  one  minute  after  the  signal  arrival  for  Event  3, 

4  1/2  seconds  of  data  were  flagged  as  invalid  and  contained  no  valid  seismic  in¬ 

formation.  A  tape  dump  of  one  of  the  7  signal  samples  not  successfully  edited 
also  showed  that  the  data  were  flagged  as  invalid  and  contained  no  valid  seismic 
information. 

The  15  signals  and  14  noise  samples  successfully  edited  appeared 
to  have  no  recording  problems  with  the  exception  that,  in  the  case  of  6  events, 
one  subarray  was  down  during  both  the  signal  and  noise  samples. 

B.  Noise  Analysis 

The  noise  analysis  is  based  on  14  five-minute  noise  samples  asso¬ 
ciated  with  the  events  of  Table  V-l  and  covering  the  period  9  March  1971  to  3  May 
1971.  Because  of  the  similarity  of  these  samples,  data  from  the  sample  asso¬ 
ciated  with  the  Greece  event  has  been  used  for  display  purposes. 

The  preliminary  noise  analysis  has  two  basic  objectives:  the  first  is 
to  examine  the  general  data  quality,  the  second  is  to  examine  noise  characteristics 
as  they  affect  array  processing. 

The  following  measurements  were  made  to  meet  these  objectives: 

•  Visual  review  of  the  data 

•  Single  -  sensor  RMS  levels 

•  Single  -  sensor  spectra  and  spectral  ratios 

•  Subarray  beam  spectra  and  spectral  ratios 

•  Array  beam  spectra 
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•  Multiple  coherence  within  subarrays 

•  Multiple  coherence  between  subarrays 

A  visual  review  of  the  noise  data  showed  that  there  were  no  spikes, 
glitches,  or  other  anomalies  in  the  samples  and  that  all  traces  appear  to  be 
"normal". 

The  single  sensor  RMS  noise  levels  for  each  of  the  132  sensors 
were  calculated  to  determine  the  distribution  of  wide-band  single  sensor  power. 
Figure  V-l  is  a  histogram  for  the  RMS  power  for  the  132  data  channels  for  the 
Greece  noise  sample  and  shows  that  the  spread  of  RMS  noise  values  is  small, 
suggesting  that  the  sensors  are  well  equalized  and  the  absolute  noise  level  is 
rather  constant  over  the  array.  Figure  V-2  is  the  corresponding  histogram  for 
a  noise  sample  taken  before  event  11  from  Southern  Iran.  The  spread  of  RMS 
values  is  similar  to  that  shown  in  Figure  V-l.  RMS  values  before  event  8  from 
the  Kurile  Islands  showed  an  even  smaller  spread. 

The  single  sensor  power  spectra  and  variations  in  these  spectra 
were  examined  next.  The  characteristics  of  the  single  sensor  and  array  noise 
field  are  summarized  in  Figure  V-3,  which  shows  the  average  single  sensor 
spectrum,  array  beam  spectrum,  and  the  spectral  ratio  of  the  array  beam  to  the 
average  single  sensor  spectrum  for  the  Greece  sample.  No  instrument  response 
corrections  have  been  made  to  these  data.  Also,  the  spectra  are  not  absolute. 
The  noise  spectra  are  very  simple,  with  a  peak  at  .  16  to  .  33  Hz  (6  to  3  second 
period)  and  a  very  rapid  decrease  in  spectral  density  with  increasing  frequency. 
Note  that  the  array  beam  achieves  ^N  noise  rejection  over  the  entire  0  to  5  Hz 
band. 

The  single  sensor  power  spectra  variations  are  summarized  in 
Figure  V-4,  which  shows  a  typical  sensor  to  average  single  sensor  spectral 
ratio,  the  maximum  positive  ratio  (sensor  3  of  subarray  10)  and  the  maximum 
negative  ratio  (sensor  2  of  subarray  11).  Figure  V-4  shows  that  the  132 
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HISTOGRAM  FOR  NORSAR  SHORT  PERIOD  NOISE  RMS  LEVELS 

GREECE  NOISE  SAMPLE 
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FIGURE  V-2 

HISTOGRAM  FOR  NORSAR  SHORT  PERIOD 
NOISE  RMS  LEVELS  BEFORE  EVENT  11 
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FIGURE  V-3 

ARRAY  BEAM  &  AVERAGE  SINGLE  SENSOR  NOISE  POWER  SPECTRA 

GREECE  N^S^  SAMPLES 


Ratio  (dB  Down) 
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single  sehsor  power  spectra,  are  very  similalr;  the  maximum  positive  or  neg¬ 
ative  variation  is  about  5  dB,  ^nd  for  most  sensors  it  is  less  than  2  dB  at  all 

i 

frequencies  i 


,  Only  one  anomaly  was  noted  in  the  single  sensor  power  spectra; 

sensors  0  and  1  of  subarray  22  had  a  spectral  line  qt  2.  8  Hz.  This  anomaly 
persisted  in  the  remaining  13  noise  samples.  i 


!  Figure  V- 5  .shows  the  array  beam  spectrum,  the  average  sub¬ 
array  beam  fepectrum,  and  the  spectral  ratio  of  the  array  beam  to  the  average 

i  1 

subarray  beam  for  the  Greece  noise  sample.  Again,  the  sinpple  spectral  shape 

.  1  !  I' 

of  the  noise  apd  the  VN~  npise  rejection  over  the  entire  band  is  evident. 


i  Variations  in  the  noise!  spectra  qf  the  subarray  beams  were  ex¬ 

amined  by  comparing  the  ratio  of  eaich  subarray  beam  to  the  average  subarray 
beapi.  In  the  region  of  significant  power  density  (0  to  2.0  Hz)  all  subarrays 
were  within  2  dB  of  the  ayerage  subarray  spectrum.  * 

l  • 

f  .  i 

Multiple  cqherences  within  a  subarray  (five  sensors  on  ring  pre¬ 
dicting  the  ceptral  sensor  of  tlje  subarray)  were  calculated  for  the  22  subarrays. 

!  !  !  . 

A  typical  multiple  Coherency  plot  (subarray  10)  is  given  in  Figure  V-6  for  vhe 

1  i  , 

Greece  sample.  Except  for  the  peak  at  3  to  6  seconds,  coherence  is  'low  at  all 
frequencies  where  there  is  significant  energy.  The  3  to  6  second  energy  is  well 
below  £he  signal  band,  and  so  can  be  removed  by.  band  pass  filtering.  Above 

2  Hz,  energy  levels  are  at  least  25  dB  below  the  spectral  peak,  so  that  there  is 

,  1  ! 

little  energy  to  remove.  1  The  moderate  coherence  here. may  possibly  be  due  to 
the  system  electronics.  In  any  e'vent  it.is  clear  that  multichannel  processing 

;  ■  • .  j  . 

would  not  be  necessary'at  the  subarray  level. 

!  .  j  J 

!  '  .  ,  ,  .  •'  ! 

Multiple  coheriencies  calculated  between  subarrays  were  very  low, 

,  i 

as  expected,  because  of  the  large  distance  involved.  Figure  V-7  shows  the 
multiple  cdherence  between  the  reference  sensor  of  subarray  1  and  the  reference 

seqsors  of,  subarrays  2,  3,  4,  6  and  8  for  the  Greece  sample;  very  low  values  were 
obtained.  !  !  y_g  i  ( 
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NOISE  POWER  DENSITY  FOR  ARRAY  AND  AVERAGE  SUBARRAY  BEAMS 

GREECE  NOISE  SAMPLE 
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FIGURE  V-6 

MULTIPLE  COHERENCE  FOR  SUBARRAY  10  NOISE  FIELD 
GREECE  NOISE  SAMPLE 
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Multiple  Coherence 
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FIGURE  V-7 

MULTIPLE  COHERENCE  BETWEEN  THE  REFERENCE 
SENSOR  OF  SUBARRAY  1  AND  THE  REFERENCE 
SENSORS  OF  SUBARRAYS  2,  3,  4,  6  AND  8 
GREECE  NOISE  SAMPLE 
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Examination  of  the  13  other  noise  samples  indicates  that,  with 
rare  exceptions  confined  to  individual  subarrays,  the  noise  behaves  very  much 
like  the  noise  preceding  event  1  from  Greece.  As  expected,  however,  power 
levels  fluctuate  from  day  to  day.  Table  V-2  gives  the  RMS  level  of  the  array 
beam  for  13  noise  samples.  The  first  three  characters  of  each  noise  sample 
name  correspond  to  the  region  of  the  event  which  followed.  The  middle  three 
digits  give  the  Julian  date,  and  the  last  two  digits  give  the  hour.  RMS  levels 
are  in  computer  counts.  A  factor  of  3  (10  dB)  variability  is  observed.  However 
there  does  not  appear  to  be  a  consistent  trend  over  the  two  month  time  period 
analyzed. 

In  summary,  this  preliminary  analysis  has  shown: 

•  The  short  period  array  appears  to  be  operating  quite  reliably 

•  The  noise  spectrum  is  quite  simple,  with  a  peak  at  3  to  6  seconds 
and  a  very  rapid  decrease  with  increasing  frequency 

•  Spectral  variability  across  the  array  is  quite  small,  less  than 
+  2  dB  for  array  beams 

•  Noise  coherence  is  low,  except  at  the  3  to  6  second  peak  (which 
is  outside  the  signal  band) 

•  VN-  noise  rejection  is  achieved  by  simple  beamforming  over  the 
entire  signal  band. 

•  Approximately  10  dB  variation  in  broadband  noise  power  was 
observed  in  1  3  noise  samples  spread  over  two  months. 

C.  Signal  Analysis 

Fifteen  signals  recorded  by  the  entire  NORSAR  array  (132  sensors) 
were  analyzed  in  order  to  study  signal  characteristics.  The  following  studies 
were  performed  for  these  events: 
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TABLE  V-2 


NOISE  RMS  LEVELS  OF  ARRAY  BEAM 


- ...  — - -"1 

Event 

Name 

RMS  Level 

RMS  Level  (dB) 

1 

GRE/068/04N 

23.00 

27.  2 

4 

KIR/082/09N 

14.91 

23.  5 

5 

KIR/082/20N 

17.  87 

25.0 

6 

KIR/083/ 20N 

15.  82 

24.0 

7 

TIB/093/07N 

12.01 

21. 6 

8 

KUR/099/15N 

11. 86 

21.  5 

9 

YUG/100/02N 

9.  51 

19.  5 

10 

KUR/101/08N 

24.  22 

27.  7 

11 

IRA/102/19N 

17.  60 

24.  9 

12 

KAZ/108/23N 

23.43 

27.  4 

13 

GRE/109/02N 

28.94 

29.  2 

14 

TAD/110/03N 

24.08 

27.  6 

15 

TIB/123/00N 

i 

11.82  | 

21.  5 
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•  Signal  similarity  investigations 

•  Signal  spectrum  measurements 

•  Noise  and  signal  spectral  comparisons 

•  Delay-time  anomalies  {deviation  from  plane  wave  propagation) 

Qualitative  estimates  of  signal  similarity  were  made  by  analysis 
of  all  sensors  for  events  1  and  2  and  of  subarray  outputs  for  all  15  events. 

Within  a  subarray,  the  waveforms  were  much  the  same  for  the  entire  length 
of  the  event.  However,  between  the  subarray  outputs,  there  was  moderate  sim¬ 
ilarity  only  for  the  first  few  cycles  of  an  event.  Significant  differences  occured 
within  a  few  seconds  after  the  signal  arrival.  As  a  result,  array  beams  formed 
from  the  22  subarray  outputs  show  a  moderate  power  drop  compared  with  the 
strongest  subarray  a  few  seconds  after  a  signal  arrival. 

Signal  amplitude  variations  between  subarray  outputs  are  striking 
(note  that  the  seismometers  appear  to  be  well  equalized).  For  the  GRE/068/04N 
event,  the  peak-to-peak/2  values  are  shown  in  Figure  V-8.  These  values 
vary  from  1999  at  subarray  12  to  589  at  subarray  6  (10.6  dB  variation).  For  this 
event,  the  southwest  quadrant  of  the  NORSAR  array  (subarrays  6,  7,  15,  16,  17,  18, 
19,  20,  and  21)  exhibits  peaks  below  1000,  while  the  northeast  section  (subarrays 
12,  13,  and  14)  has  peaks  above  1500.  For  the  URA/082/06N  event,  the  corres¬ 
ponding  zero -to -peak  values  are  shown  in  Figure  V-9.  These  values  vary 
from  7830  at  subarray  11  to  5438  at  subarray  7  (3.  2  dB  variation).  Subarray  12 
still  had  the  highest  signal-to-noise  ratio  for  this  event.  In  the  old  Plan  D  events 
from  the  Kazakh  test  site,  the  signal  amplitude  variations  consistantly  exceeded 
the  variations  of  the  GRE/068/04N  event.  In  most  of  the  15  NDPC  events  pro¬ 
cessed,  subarrays  11  or  12  were  the  strongest.  However,  for  the  two  events 
from  the  Kurile  Islands  region,  subarray  3  was  the  strongest  subarray.  Since 
subarray  3  is  located  on  the  inner  ring  of  the  NORSAR  array,  it  is  clear  that  the 
location  of  the  strongest  subarray  is  no  simple  function  of  event  epicenter. 
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The  signal  spectral  characteristics  for  the  GRE/068/04N  event 
were  examined  using  the  subarray  12  beam,  which  had  the  largest  signal-to  - 
noise  ratio.  The  subarray  12  signal  plus  noise  spectrum  for  the  128- second 
gate  is  shown  in  Figure  V-10;  the  noise  spectrum  for  320-second  gate  just 
preceding  the  signal  also  is  plotted  for  reference.  The  signal  energy  is  con¬ 
centrated  in  a  band  from  0.  5  to  1.  5  Hz  with  a  sharp  drop  in  power  density  on 
either  side  of  this  band.  The  signal -to -noise  ratio  is  best  between  1.0  and  1.  5 
Hz,  suggesting  that  signal  detection  may  best  be  accomplished  at  NORSAR  at 
relatively  high  frequencies. 

Figure  V-ll  shows  the  effect  of  plane-wave  beamforming  on 
signal  strength  for  the  GRE/068/04N  event.  The  average  single  sensor  spectrum, 
the  average  subarray  beam  spectrum  (\xsing  plane  wave  delays)  and  the  array 
beam  (also  using  plane  wave  delays)  are  shown;  a  128-second  signal  gate  was 
used  to  obtain  the  spectral  estimates.  The  corresponding  noise  spectra  com¬ 
puted  from  a  320-second  signal  gate  just  preceding  the  signal  were  shown  pre¬ 
viously;  note  that  the  signal-to -noise  ratio  for  this  event  is  only  fair. 

At  the  1.0  Hz  signal  peak,  signal  loss  resulting  from  subarray 
beamforming  is  2  dB;  this  is  slightly  higher  than  expected.  However,  the  low 
signal-to -noise  ratio  at  several  of  the  subarrays  could  account  for  part  of  the 
loss.  It  appears  that  plane  wave  delays  are  adequate  for  forming  subarray 
beams,  as  is  the  case  at  LASA  (Harley,  1967). 

Plane  wave  delays  clearly  are  inadequate  for  array  beamforming. 
Signal  loss  at  1.0  Hz  is  9  dB  (relative  to  the  average  subarray  beam)  and  the 
array  beam  signal -to -noise  ratio  is  not  much  better  than  that  for  the  average 
subarray  beam  (and  in  fact  significantly  worse  than  that  for  certain  high- amplitude 
subarrays).  Thus  time  delay  anomaly  corrections  are  essential  for  success¬ 
ful  array  beamforming;  this  is  discussed  in  more  detail  in  the  next  section. 
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FIGURE  V-10 

SUBARRAY  12  SIGNAL  SPECTRUM  FOR  GRE/068/04N 
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Figure  V-12  shows  a  spectrum  for  the  first  6.4  seconds  after  the 
arrival  of  URA/082/06N  and  represents  the  power  density  for  the  output  of  a 
diversity-stack  beam.  The  bodywave  magnitude  for  this  event  is  5.  6,  and  the 
signal -to -noise  ratio  is  very  large.  In  comparison  with  the  other  14  events,  the 
spectrum  is  relatively  uniform  from  0.  5  Hz  to  beyond  2,  0  Hz,  where  it  begins  its 
final  falloff.  A  comparison  with  typical  noise  spectra  shows  that  the  highest 
signal -to -noise  ratio  for  this  event  occurs  around  2.0  Hz. 

Figure  V-13  shows  a  spectrum  for  the  first  12.8  seconds  after 
the  arrival  of  KIR/082/20N  on  a  diversity-stack  beam.  The  event  bodywave 
magnitude  is  6.0,  and  the  signal -to -noise  ratio  is  very  large.  The  spectral 
values  below  0.  5  Hz  represent  power  density  after  application  of  a  0.  5-5.  0  Hz 
bandpass  filter.  This  spectrum  is  typical  of  all  events  except  the  Ural  Mountains 
event  just  discussed.  In  the  Kirgiz  event,  the  spectrum  tends  to  fall  off  beginning 
at  1  Hz.  The  highest  signal -to -noise  ratio  for  this  event  occurs  around  1.  5  Hz. 

In  summary,  this  preliminary  signal  analysis  has  shown: 

•  Signal  similarity  appears  to  be  fairly  good  within  a  subarray,  and 
reasonably  good  in  the  first  few  cycles  of  the  signal  between  sub¬ 
arrays. 

•  Signal  amplitude  variations  across  the  array  are  significant 

•  For  all  but  the  Ural  Mountains  event,  signal -to -noise  ratios  were 
highest  between  1.0  and  1.5  Hz.  For  that  event,  the  signal-to- 
noise  ratio  was  highest  at  about  2.0  Hz. 

•  Theoretical  plane  wave  delays  appear  to  be  adequate  for  subarray 
beamforming,  but  time  delay  anomaly  corrections  are  required  for 
array  beamforming. 
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DIVERSITY -STACK  SIGNAL  SPECTRUM  FOR  6.4  SECONDS 
AFTER  ARRIVAL  OF  URA/082/06N  ! 
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D.  Array  Beamforming  Results 

As  mentioned  in  the  previous  subsection,  there  were  two  principal 
Pr  bJ  ems  encountered  in  forming  beams  to  enhance  signals  propagating  across 
the  NORSAR  array.  First,  the  wavefront  arrival  time  at  the  individual  sub¬ 
arrays  were  not  consistent  with  plane  wave  propagation  at  the  expected  velocity 
and  azimuth.  Second,  large  variations  in  the  received  signal  amplitude  occurred 
between  subarrays. 

In  an  attempt  to  remove  the  inter- subarray  delay  anomalies,  cross¬ 
correlation  functions  were  computed  between  a  reference  subarray  beam  and  the 
remaining  subarray  beams  in  order  to  obtain  an  adusted  set  of  delays.  The  ref¬ 
erence  subarray  generally  was  chosen  to  be  the  one  with  the  highest  signal-to- 
noise  ratio.  This  method  worked  well  when  the  subarrays  had  high  signal-to- 
noise  ratios,  but  for  the  subarrays  with  the  weakest  signals,  results  were  un¬ 
satisfactory.  For  these  subarrays,  delays  were  obtained  by  a  combination  of 
visual  inspection,  examination  of  crosscorrelation  functions,  and  comparison  of 
delay  anomalies  from  other  nearby  events. 

•After  removing  the  delay  anomalies,  the  subarrays  were  weighted 
to  give  a  diversity- stack  beam.  In  this  technique,  each  subarray  is  weighted 
proportionate  to  the  subarray  signal-amplitude  to  noise-power  ratio  and  added  to 
the  other  weighted  and  shifted  subarray  traces.  A  detailed  study  of  a  five-minute 
noise  sample  just  before  GRE/068/04N  event  showed  that  it  was  necessary  to 
use  noise  gates  up  to  60  seconds  in  length  before  the  power  variations  with  time 
were  reduced  to  the  same  level  as  the  variations  in  power  between  the  subarray 
outputs  over  the  entire  five-minute  noise  sample.  Based  on  the  noise  analysis 
results  shown  previously  the  assumption  of  equal  noise  power  on  all  subarray 
beams  is  quite  reasonable.  All  subarray  outputs  were  therefore  assumed  to  have 
the  same  noise  power.  The  signal  amplitude  for  each  subarray  was  computed 
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using  the  RMS  level  over  a  6.4-second  or  12.  8-second  window  beginning  with 
the  signal  arrival. 

Beams  were  formed  in  this  way  for  six  of  the  15  signals  success¬ 
fully  edited.  Results  for  these  events  are  shown  in  Table  V-3  in  terms  of  the 
resulting  signal-to-noise  ratios  for  the  reference  subarray,  the  plane-wave 
beam,  the  adjusted- delay  beam,  and  the  diversity-stack  beam.  Noise  power 
was  computed  over  a  70-second  gate  preceding  the  signal.  Signal  power  was 
computed  over  a  6.  4-second  or  12.8-second  window  beginning  with  the  signal 
arrival.  The  reference  subarray  used  wa3  always  the  subarray  with  the  highest 
signal-to-noise  ratio.  A  common  set  of  adjusted  delays  was  applied  to  the  three 
Kirgiz-Sinkiang  border  events. 

Improvement  in  the  diver sity- stack  beams  ranged  from  2.  9  to  8.  1  dB 
relative  to  the  reference  subarray,  from  1.  7  to  9.  3  dB  relative  to  the  plane-wave 
beam,  and  from  -0.  1  to  2.  3  dB  relative  to  the  adjusted-delay  beam.  Most  of  the 
improvement  is  obtained  from  application  of  the  adjusted  delays.  In  only  one 
case  did  the  improvement  over  the  adjusted-delay  beam  exceed  0.  7  dB.  The  use- 
fullness  of  the  diversity-stack  technique  will  require  additional  analysis.  In  part¬ 
icular,  a  small  ensemble  of  events  from  Eastern  Kazakh  will  be  used  for  eval¬ 
uation.  The  old  "Plan  D"  data  indicated  that  signal  amplitude  variations  were 
largest  for  this  location,  and  hence  the  diversity-stack  technique  would  be  ex¬ 
pected  to  give  the  largest  improvement  over  conventional  beamsteering,,. 

Plane-wave  and  adjusted  delays  for  Eastern  Kazakh,  Alma-Ata, 
and  the  GRE/068/04N  event  are  shown  in  Table  V-4.  The  Eastern  Kazakh 
and  Alma-Ata  delays  were  obtained  from  "Plan  D"  data  recorded  in  1970.  The 
Greece  event  had  only  a  moderate  Signal-to-noise  ratio,  so  that  the  tabulated 
delay  anomalies  are  less  certain  than  for  the  other  events.  Delay  anomalies  of 
0.  5  seconds  were  typical;  note  that  an  anomaly  of  this  size  is  enough  to  cause 
half-cycle  errors  in  the  time  alignment  of  1  He  signals  (and  hence  cancellation). 
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TABLE  V-3 


BEAM  SIGNAL -TO -NOISE  RATIOS  FOR  SIX  EVENTS 


Event 

Bandpass 

(Hz) 

5 

mb 

T - 

Signal-to 

-Noise  Ratio  (dB) 

Reference 

Subarray 

Planc- 

i 

Wave 

Beam 

I  Adjusted- 
I  Delay 
Beam 

Diversity 
Stack  1 
Beam  i 

URA/082/06N 

0.05-5.0 

5.6 

38.  4 

39.  3 

43.  1 

- 

43.  6 

,KIR/082/09N 

0.  5-5.0 

5.7 

38.  2 

37.0 

40.8 

1 

4i.  t  ; 

KIR/082/20N 

0.  5-5.0 

6.0 

32.4 

36.  7 

38.  5 

i 

38. 4  t 

KIR/083/ 20N 

0.  5-5.0 

5.  3 

30.  1 

30.  2 

35.8 

36.  5 

TIB/093/07N 

0.  5-5.0 

5.  1 

6.  3 

6.  1 

14.  1 

( 

14.4 

YUG/100/02N 

0.  5-5.0 

4.  6 

18.6 

13.  1 

20.0 

22.4 
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using  the  RMS  level  over  a  6.4-second  or  12.8-second  window  beginning  with 
the  signal  arrival. 

Beams  were  formed  in  this  way  for  six  of  the  15  signals  success¬ 
fully  edited.  Results  for  these  events  are  shown  in  Table  V-3  in  terms  of  the 
resulting  signal-to-noise  ratios  for  the  reference  subarray,  the  plane-wave 
beam,  the  adjusted- delay  beam,  and  the  diversity-stack  beam.  Noise  power 
was  computed  over  a  70-second  gate  preceding  the  signal.  Signal  power  was 
computed  over  a  6.  4-second  or  12.  8- second  window  beginning  with  the  signal 
arrival.  The  reference  subarray  used  was  always  the  subarray  with  the  highest 
signal-to-noise  ratio.  A  common  set  of  adjusted  delays  was  applied  to  the  three 
Kirgiz-Sinkiang  border  events. 

Improvement  in  the  diversity-stack  beams  ranged  from  2.  9  to  8.  1  dB 
relative  to  the  reference  subarray,  from  1. 7  to  9.  3  dB  relative  to  the  plane-wave 
beam,  and  from  -0.  1  to  2.  3  dB  relative  to  the  adjusted-delay  beam.  Most  of  the 
improvement  is  obtained  from  application  of  the  adjusted  delays.  In  only  one 
case  did  the  improvement  over  the  adjusted-delay  beam  exceed  0.  7  dB.  The  use- 
fullness  of  the  diver sity- stack  technique  will  require  additional  analysis.  In  part¬ 
icular,  a  small  ensemble  of  events  from  Eastern  Kazakh  will  be  used  for  eval¬ 
uation.  The  old  '  Plan  D"  data  indicated  that  signal  amplitude  variations  were 
largest  for  this  location,  and  hence  the  diversity-stack  technique  would  be  ex¬ 
pected  to  give  the  largest  improvement  over  conventional  boamsteering. 

Plane-wave  and  adjusted  delays  for  Eastern  Kazakh,  Alma-Ata, 
and  the  GRE/068/04N  event  arc  shown  in  Table  V-4.  The  Eastern  Kazakh 
and  Alma-Ata  delays  were  obtained  from  "Plan  D"  data  recorded  in  1970.  The 
Greece  event  had  only  a  moderate  signal-to-noise  ratio,  so  that  the  tabulated 
delay  anomalies  arc  less  certain  than  for  the  other  events.  Delay  anomalies  of 
0.  >  seconds  were  typical;  note  that  an  anomaly  of  this  size  is  enough  to  cauto 
half-cycle  errors  in  the  time  alignment  of  1  Hz  signals  (and  hence  cancellation). 
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The  close  agreement  between  the  delay  anomalies  for  the  Kazakh 
area  and  the  Alma-Ata  event  (separated  by  7.  5  degrees  in  latitude  and  0.  6  degrees 
in  longitude)  indicates  that  a  single  set  of  time  delay  corrections  may  be  useful 
for  a  relatively  large  region.  Delay  anomalies  obtained  for  the  events  taken 
from  the  NORWAY  NDPC  tapes  tend  to  confirm  the  fact  that  the  delay  anomalies 
change  relatively  slowly  as  a  function  of  event  epicenter.  The  delay  anomalies 
for  the  GRE/068/04N  event,  for  example,  were  useful  in  deciding  the  delay 
anomalies  for  the  GRE/109/02N  and  YUG/100/02N  events. 

Figure  V-14  shows  the  URA/082/06N  event  and  illustrates  the 
improvement  in  signal -to-noise  ratio  resulting  from  time  delay  corrections  and 
the  diversity-stack  technique.  The  top  trace  is  the  plane-wave  beam,  the  second 
trace  is  the  adjusted-delay  beam,  the  third  trace  is  the  diversity-stack  beam, 
and  the  fourth  trace  is  subarray  12,  the  reference  subarray.  All  traces  have 
been  filtered  to  pass  the  band  0.05  to  5.0  Hz.  The  traces  have  been  scaled  so 
that,  if  all  subarray  outputs  sum  in  phase,  the  resulting  beams  would  have  iden¬ 
tical  ampltudcs. 

Due  to  the  fact  that  the  delays  used  in  the  plane-wave  beam  arc 
inaccurate,  waveform  distortion  and  signal  attenuation  arc  evident  in  the  top 
wave  beam  relative  to  the  reference  subarray  shown  in  the  bottom  trace.  Use  of 
the  adjusted  delays,  however,  produced  a  significant  4  dB  improvement  in  the 
adjusted-delay  and  diversity-stack  beams  relative  to  the  plane-wave  beam.  Be¬ 
cause  of  the  lack  of  signal  similarity  after  the  first  seconds  of  the  signal,  the 
two  middle  traces  show  a  marked  power  drop  as  compared  with  the  reference 
subarray  at  the  bottom. 

E.  Future  Plans 

During  the  third  quarter  of  this  contract,  routine  processing  of 
events  will  continue.  In  addition,  events  processed  will  be  run  through  the  Texas 
Instruments  discrimination  subroutine. 
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VI. 


LONG  PERIOD  EXPERIMENT 


A.  Introduction 

This  section  describes  the  work  performed  during  the  past  quarter 
on  the  Long  Period  Experiment  (LPE).  The  overall  objectives  of  our  LPE  eval¬ 
uation  are  to  determine  the  single  station  and  multistation  detection  threshold  and 
to  determine  the  multistation  discrimination  capability.  To  date  the  evaluation  has 
been  limited  by  the  quality  and  quantity  of  the  data  received  at  SAAC.  This  limit¬ 
ation  was  a  result  of  problems  with  unreliable  station  power  supplies,  component 
polarity  reversals,  and  unequalized  seismometer  responses.  However,  Lamont 
recently  has  sent  several  crews  to  the  stations  to  correct  these  problems,  and  it  is 
expected  that  better  data  should  be  available  in  the  near  future  from  the  first  five 
LPE  stations.  (Note  that  data  recently  was  received  from  .hree  of  the  stations 
after  the  modification  period).  Table  VI-1  lists  all  data  received  at  SAAC. 

Because  of  the  problems  outlined  above,  only  limited  amounts  of 
data  were  processed  during  the  quarter.  The  times  processed  cover  the  periods 
January  24,  1971  to  February  6,  1971  for  stations  (sites)  1,2,  3,5,7  (Table  VI-2) 
and  February  9,  1971  to  March  15,  1971  for  station  3.  This  section  describes 
the  results  obtained  from  analyzing  these  time  periods.  Included  arc  an  analysis 
of  the  spectral  content  and  time  variability  of  the  noise  field  at  each  station,  a 
very  preliminary  discussion  of  the  detection  capability  of  each  station  and  the 
network,  and  the  behavior  of  the  Ms-mjj  classification  parameter. 

B.  Noise  Analysis 

The  noise  analysis  discussed  in  this  section  covers  data  from  the 
time  period  January  24,  1971  to  February  6,  1971.  Data  is  presented  for  the 
five  stations  in  Table  VI-2.  It  should- be  noted  that  data  discussed  in  this  section 
was  not  of  the  highest  quality,  and  that  it  does  not  represent  the  current  record¬ 
ing  capability  of  the  stations.  Therefore,  these  results  may  not  be  representative 
of  the  current  noise  field. 
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TABLE  VI- 1 

LPE  DATA  RECEIVED  AT  SAAC  FROM  LAMONT  GEOLOGICAL 

OBSERVATORY 

All  days  arc  Julian  (Only  data  from  1971  is  shown) 


AUSTRALIA 


THAILAND 


ALASKA 


ISRAEL 

OGDENSBURG 

1-4 

7-55 

1-37 

42-260 

1-9 

21-109 

220-236 


TABLE  VI-2 

STATIONS  OF  THE  LPE  NETWORK 


Station 

Code 

Coordinates 

Elevation 

(ft) _ 

...  - 

Latitude 

(deg: min:  sec) 

Longitude 
(deg:  min:  see) 

Australia 

1 

20:05:  18  S 

146:15:16  E 

1168 

Thailand 

2 

18:47:24  N 

98:  57:  37  E 

1352 

Alaska 

3 

64:53:58  N 

148:00:20  W 

1000 

Israel 

5 

29:19:48  N 

34:  32:  12  E 

655 

Ogdensbi  rg 

7 

41:04:00  N 

74:37:00  W 

-1221 

1. 


Spectral  Content 


To  study  the  spectral  content  of  the  noise,  a  one-hour  noise  sample 
from  each  day  of  the  time  period  was  edited.  However,  since  visual  records 
of  the  data  were  not  available  when  the  sample  time  periods  were  chosen,  several 
of  the  samples  contain  unwanted  signals.  The  time  periods  containing  signals  or 
questionable  data  are  not  included  in  the  noise  analysis.  Also,  as  will  be  seen 
in  the  figures  of  this  and  the  signal  analysis  section,  none  of  the  stations  was 
operative  throughout  the  entire  14  days  of  interest. 

Figure  Vl-1  shows  noise  power  density  spectra  for  two  days. 

These  spectra  have  been  corrected  to  millimicrons  of  ground  motion  using  cal¬ 
ibration  information  furnished  by  Lamont  Observatory.  (The  corrections  arc 
relative  to  l(mjj  )^/H z  at  0.025  Hz).  In  the  figure  only  the  vertical  components 
have  been  presented  to  show  the  relative  noise  power  at  each  of  the  stations.  This 
figure  is  typical  of  the  spectra  for  the  time  period  studied  and  exhibits  certain 
properties  of  the  five  noise  fields  which  arc  discussed  below. 

The  noise  fields  for  sites  2  and  3  arc  dominated  by  a  large  peak  at 
the  very  long  periods.  This  peak  is  present  on  all  of  the  spectra  analyzed  for 
the  two  sites.  Site  1  also  contains  a  long  period  peak  for  all  but  one  of  the  samples; 
however,  it  is  not  as  dominant  as  the  peak  for  sites  2  and  3.  Site  7  contains  a  long 
period  peak  on  less  than  half  of  the  samples  analyzed;  note  that  in  the  two  spectra 
shown,  site  7  has  an  extremely  large  peak  in  the  bottom  spectra,  but  the  lop 
spectra  is  relatively  flat  out  to  0.07  Hz.  Also,  site  5  docs  not  contain  a  long 
period  peak  in  any  of  the  samples  analyzed;  however,  the  data  for  this  site  were 
limited  to  the  time  period  Janaury  24,  1971  to  January  29,  1971. 

For  typical  long  period  noise  fields  the  systems  were  designed  so 
that  the  oise  spectra,  as  seen  through  the  systems,  will  be  flat  in  the  20-50 
second  range.  This  property  was  observed  in  some  of  the  spectra  for  sites  5  and7. 
(Figure  VI-2) 
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'  1  '  , 

<  2.  Time  Variability 

i  ;  i 

,  The  time  variability  of  the  noise  field  at  each  of  the  sites,  over  the 

January  24,  1971  to  February  6;  1971  time  period  was  studied  by  computing  the 
RMS  level  in  the  20-40  second  band.  After  computation,  the  levels  were  correct- 

I  «  .  f  * 

ed  to  absolute  ground  motion  using  the  Lamont  calibration  data.  It  should  be 
noted  that  the  20-40  second  responses  are  reasonably  close  between  the  vertical 
and  horizontal  (components  for  sites  1,2,  and  5.  However,  the  vertical  response 
for  site  3  (Alaska)  is  about  twice  that  of  the  North-South  component.  Thus  a 

1  .  i 

direct  component  comparison  can  be  made  for  all  sites  except  3  (Note:  no  responses 
were  available  for  site  7).  RMS  calculations  are  shown  in  Figure  VI-3.  From 

:  '  1  I 

this  figure  the  following  conclusions  have  been  made: 

•  i 

Site  1  -  very  little  can  be  said  due  to  the  limited  amount  of  data  but 
the  East-West  horizontal  appears  to  be  the  noisiest.  The 
difference  between  the  maximum  and  minimum  component  RMS 
levels  arc:  vertical  =  4.4,  North-South  =2.7,  and  East-Wcst= 

1  6.0  ' 

Site  2  -  an  increasing  trend  in  noise  level  with  time  is  seen  for  all 

*  i 

three  components.  The  difference  between  the  maximum  and 

*  J  , 

minimum  component  RMS  levels  arc:  verticals  4.  5,  North- 
South=6.  6,  and  East-West=  3.  8.  Also,  the  two  horizontal 
components  arc  much  higher  than  the  vertical  component.  This 

,  » 

difference  also  is  visible  in  the  spectra  and  is  most  pronounc- 

I 

i  cd  at  the  longer  periods.  i 

Site  3  -  as  opposed  to  site  2,  the  vertical  component  is  the  noisiest  for 
site  3.  A  possible-explanation  for  the  noisier  vertical  is  that 
the  conversion  to  ground  motion  was  done  using  the  calibration 
information  at  40  sccorids.  Because  of jthc  difference  in  am¬ 
plitude  response  between  the  vertical  and  horizontals  for  Alaska, 

i 

I 

I  ' 

‘  '  •  I  |  I 
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Tlm^ 

FIGURE  VI- 3 
NOISE  RMS  LEVELS 


the  vertical  would  l  ave  a  higher  RMS  level  than  the  horizontals 
in  the  20-40  second  band  if  the  absolute  noiuc  spectra  were 
identical.  However,  the  noise  levels  do  tend  to  remain  fairly 
constant  through  most  of  the  time  period  (Note:  the  North- 
South  component  was  not  operating  correctly  for  samples  1-4). 

The  difference  between  the  maximum  and  minimum  component 
RMS  levels  are:  verticals  1.  7,  North-South=  3.  1,  and  East- West=  0.  9 
the  noise  level  for  the  vertical  and  North-South  components 
remain  fairly  level.  However,  observing  the  large  difference 
between  the  two  horizontal  RMS  levels,  it  appears  that  the  am¬ 
plitude  scale  factor  may  be  suspect  for  the  East-West  compon¬ 
ent.  The  difference  between  the  maximum  and  minimum  com¬ 
ponent  RMS  levels  arc:  verticals  1.  2,  North-South- 2.  5,  and 
East-Wcst=  3.  7. 

the  vertical  and  East-West  components  appear  to  remain  fairly 
constant  throughout  the  time  period.  The  difference  between 
the  maximum  and  minimum  component  RMS  levels  arc: 
verticals  1.  6,  and  East- Wests 0.  9.  Again  the  vertical  com¬ 
ponent  is  noisier.  Since  no  frequency  response  was  available 
for  site  7,  we  could  not  check  to  sec  if  the  noisier  vertical 
component  was  due  to  the  instruments  responses  as  for  site  3. 

(Note  -  the  North-South  component  appeared  to  have  been  re¬ 
placed  with  another  vertical  instrument,  since  it  had  essen¬ 
tially  the  same  data  points  as  the  vertical  component. 


In  general,  except  for  site  2,  the  noisx:  field  RMS  levels  tend  to  remain  fairly 
constant  during  the  time  period.  Also  the  RMS  levels  for  sites  2  and  3  vertical 
components,  site  5  vertical  and  North-South  components,  and  site  7  East-West 
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component  all  are  similar.  Further,  it  appears  that  the  horizontal  components 
are  noisier  than  the  vertical  components.  The  LPE  RMS  levels  are  the  same 
order  of  magnitude  as  those  observed  at  ALPA  and  NORSAR. 

C.  Signal  Processing 

1.  Event  Detection 

This  section  discusses  the  detection  capability  of  the  five  stations 
and  of  the  network  for  the  19  events  (Table  VI-3)  recorded  during  the  period 
January  24,  1971  to  February  6,  1971.  This  event  ensemble  includes  all  reported 
events  from  the  Asian  continent  during  the  time  period  (both  USCGS  and  LASA), 
except  those  which  were  buried  in  la~ge  interfering  events  at  all  sites.  Also,  any 
site  which  was  down  or  the  data  was  bad  is  not  included  in  the  detection  or  not- 
detected  plots.  All  results  are  from  the  vertical  component  unless  otherwise 
specified,  and  all  results  are  based  on  data  which  has  only  been  bandpass  filtered 
(no  matched  filtering  or  2-component  processing  was  performed). 

Figure  VI-4  is  a  detection  plot  for  m^  versus  delta  -  the  following 
observations  are  made  from  this  figure: 

Site  1  -  all  events  with  m^  >  4.  8  were  detected 

Site  2  -  all  events  with  delta  <  50°  and  all  events  with  mjj)  4.  8  were 

detected 

Site  3  -  all  events  with  delta  <  60°  and  all  events  with  m^  y  4.  8  were 
detected 

Site  5  -  all  events  with  delta  (  80°  and  all  events  with  m^  4.  8  were 
detected 

Site  7  -  all  events  with  delta  (  75°  and  all  events  with  m^  4.  8  were 
detected 

The  dotted  lines  on  Figure  VI-4  separate  the  plots  into  zones  of  com¬ 
plete  detection  (for  the  data  sample  analyzed)  and  partial  detection;  note  that  several 
events  with  m^<[4,  8  and  delta ^>70°  were  detected  at  sites  3,  5,  and  7. 


VI-10 


h 

Z 


\Wpq«22;  \  'ST\\  2  \  ffl 


r^t^t^r^r^t'-oooot^ot'-o' 


<s  S  «  «\  z\5 

rt  * 

h 

to  oovomvomoN^ 

•-*  ^  OOOOOOOOOOhQOhOO 


nJ 

X 

CO 

10  rt 


2\  2 \W  2  \  «  2  W 


h  in  O'  in  co  oo  in  o  •— i 

t'-rot'-iMf'-t'-sO  oo  Nn- 


co 

vO 

c 

CO 

Eh 

3 

1 

HH 

2 

i'j 

i—i 

a 

u 

> 

W 

0 

an 

6  c 
O  o 


C*5  <u 
(J  T3 


tS*  \NfflfPcqzz\  . ^ \  >  *  \ 


J^22 


lOlOiOO^lOvOcOO^  OO  (\J  vO  ^  00  0s  I— I 

N  lololololololo^  ro  vD  rvj  ^  i.n  in 


i  w  (p  «  n  2  2  2 
Q 

JK  OOOinOvOOOCT' 

^  t^i^-r^r-r^^vOo 


w 


W  3  ci 


CT'O'Or^r^oooOvOvOvomoo 


wwwwwwwwwwww 

fOHifl'f^NOOKlHNl'l 


if  vO  in  5  m  m  in  o 

.  .  .  <2  .  .  w  .  ” 

J  J 


vOvOvO(V]\OOsTt<f«-POOsOLn 

inininvOinpoin^Nor'pooo 


NNHOvNCOONors>LO(MrO 


w 

W 

W 

W 

W 

w 

w 

00 

00 

O' 

0 

CO 

O' 

H 

00 

ifi 

• 

O' 

00’ 

• 

CO 

• 

O' 

O' 

'O 

O' 

t*- 

vO 

2  2  2 

00  O  o 


CT'CT'CT'roC'roi'ininor^no 

^^’tin^^^minrom'f 


u  § 

w  I  2 

tH  U  -  * 

w  «  Q 

Q  <u  o  o 

Q  2  2 


^  2  I 


PO  PO  T}< 

n  in  ^ 


W  u 
2 


'O'tvOpninHinNciHMn  r^iDvO 

HTfm^Npnino^^N^  ino^ 

h  aj  t»  h  n  #  n-  h  tn  in  ni  h  <j\  [£.’ 

in  in  (M 

pnpnpnocoi'iO'in  hoh^  p-  oi  in 

hhhqhNhhNNOh  OOO 


NnllMNNNNNNmHH  pj  rf 

- - -  ^  V  Vs 

hhhhhhhhhhNN  POIMPJ 


afflEEf 


5  bj1?""  e" 


MS™ 


::k::i": 

ssiss: 

mm 

tftssiii 


mcs 


aawroroaat 

e^ishe:: 


n*n 


EVENT  DETECTION  (nw.  v*«  Dalla}  FOR  INDIVIDUAL  SITES 


The  data  of  Figure  VI-4  were  replotted  with  Ms  vs  delta,  where 
Ms  was  calculated  by: 


Ms  -  log  “  + 

.  92  +  log  A 

A<  25° 

w  ,  A 

Ms  =  log  —  + 

1 . 66  log  A 

A^25° 

(Evernden,  1971) 

Note  that  two  of  the  events,  m^=  4.  5  and  4.6,  were  not  detected  at  any  stations 
and  thus  no  Ms  value  for  then  could  be  calculated.  For  a  given  event  Mg  was 
averaged  over  all  detecting  stations  and  this  average  value  was  used  to  plot  Ms 
versus  delta  for  each  site.  The  Ms  versus  delta  plot  for  the  individual  sites  is 
shown  in  Figure  VI-5.  Observations  from  this  figure  are: 

•  The  events  tend  to  separate  into  detected  and  not-detected  classes, 
although  there  is  some  overlap  (sites  2,  3,  and  7). 

•  There  is  a  tendency  for  the  "  detection  threshold"  to  increase  with 
distance  (as  would  be  expected)  but  the  sparse  data  are  not 
sufficient  to  specify  a  "detection  threshold  line" 

•  It  appears  that  at  60°  most  stations  have  a  detection  threshold  be¬ 
tween  Mg=  3.  0  and  Mg  =  3.  4.  The  small  range  of  this  detection 
level  among  stations  agrees  with  the  small  variance  in  RMS  noise 
levels  for  their  vertical  components.  (Note:  there  is  not  enough 
nois'e  data  to  make  this  comparison  for  site  1). 

Figure  VI-6  is  a  composite  of  Figures  VI-4  and  VI-5.  The  top 
half  of  the  figure  shows  m^  versus  delta  for  all  stations.  As  in  Figure  VI-4, 
scattering  is  evident,  however,  for  the  events  analyzed  the  following  observations 
can  be  made  for  the  five -station  network: 

•  all  events  with  m^  )  4.  8  were  detected  at  all  operating  stations 

•  all  events  were  detected  at  all  operating  stations  having  delta  <  50° 

The  bottom  half  of  the  figure  is  a  plot  of  delta  versus  Mg  for  all 
stations.  This  plot  shows  the  expected  increase  in  detection  threshold  with  dis- 
tistance.  The  broken  lines  show  the  probable  range  in  detection  threshold  as  a 
function  of  distance*  because  all  stations  appear  to  have  similar  sensitivity  this 
curve  is  a  very  crude  estimate  of  an  individual  station  detection  capability. 
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FIGURE  VI -6 


NETWORK  DETECTION 
VI- 15 


It  is  important  to  remember  that  these  graphs  are  not  intended  to 
present  a  definitive  network  or  individual  station  detection  threshold,  but  are 
presented  only  to  show  the  results  obtained  from  analysis  of  a  very  limited 
number  of  events. 

2.  Behavior  of  Standard  Discriminant 

To  study  the  behavior  of  a  standard  discriminant  for  the  LPE  net¬ 
work,  an  Ms  versus  mj,  plot  for  the  individual  sites  was  prepared.  (Figure  VI-7). 
This  figure  shows  a  large  amount  of  scattering  in  the  Ms  values  for  any  given 
event  (as  much  as  1.2  units  for  the  5.  3  mb  event).  Note  that  each  event  is  coded 
with  the  same  symbol  for  all  sites.  The  Gutenburg  Ms-mb  relationship  for 
earthquakes  (Richter,  1958)  fils  the  observed  data  reasonably  well.  All  of  the 
events  are  separated  from  the  broken  line,  which  was  obtained  from  Capon's 
data  (1967)  for  a  series  of  Asian  explosions. 

Figure  VI-8  is  a  repeat  of  Figure  VI-7  using  an  average  Ms  for 
each  event.  It  should  be  noted  that  the  Ms=  2.  6  event  was  calculated  from  one 
site  and  that  is  was  calculated  from  a  probable  LR  phase;  an  S  phase  from  the 
event  was  clearly  detectable  on  the  LP  vertical  component,  while  the  LR  wave 
was  not  larger  than  the  largest  noise  peaks.  However,  again  we  see  a  fairly 
good  fit  to  the  Gutenburg  curve  and  generally  good  separation  from  Capon's  curve. 

D.  Future  Plans 

In  the  near  future  we  expect  to  begin  receiving  data  from  all  five 
of  the  stations  in  Table  VI-1  plus  the  station  in  Spain.  When  these  data  are  received 
we  plan  to  enter  a  routine  processing  mode  (20-40  events  per  month).  However, 
until  such  data  are  available,  we  plan  to  begin  processing  the  data  for  August, 1971 
(after  the  station  modifications)  for  Australia,  Alaska,  and  Ogdensburg.  Current 
data  at  SAAC  for  the  period  is  for  days  220-236  (August  8-24). 
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